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ABSTRACT 

We present the results of a photometric redshift analysis designed to identify redshift z > 6 
galaxies from the near-IR HST imaging in three deep fields (HUDF, HUDF09-2 & ERS) 
covering a total area of 45 sq. arcmin. By adopting a rigorous set of criteria for rejecting 
low-redshift interlopers, and by employing a deconfusion technique to allow the available 
ultra-deep IRAC imaging to be included in the candidate selection process, we have derived a 
robust sample of 70 Lyman-break galaxies (LBGs) spanning the redshift range 6.0 < z < 8.7. 
Based on our final sample we investigate the distribution of UV spectral slopes (f\ oc A^), 
finding a variance-weighted mean value of (/3) = —2.05 ± 0.09 which, contrary to some 
previous results, is not significantly bluer than displayed by lower-redshift starburst galaxies. 
We confirm the correlation between UV luminosity and stellar mass reported elsewhere, but 
based on fitting galaxy templates featuring a range of star-formation histories, metallicities 
and reddening we find that, at z > 6, the range in mass-to-light ratio (M*/ Luv) at a given UV 
luminosity could span a factor of ~ 50. Focusing on a sub-sample of twenty-one candidates 
with IRAC detections at 3.6/im we find that L* LBGs at z ~ 6.5 have a median stellar mass 
of M* = (2.1 ± 1.1) x 1O 9 M0 (Chabrier IMF) and a median specific star-formation rate 
(sSFR) of 1.9 ± 0.8 Gyr" 1 . Using the same sub-sample we have investigated the influence 
of nebular continuum and line emission, finding that for the majority of candidates (16 out of 
21) the best-fitting stellar masses are reduced by less than a factor of 2.5. However, galaxy 
template fits exploring a plausible range of star-formation histories and metallicities provide 
no compelling evidence of a clear connection between star-formation rate and stellar mass at 
these redshifts. Finally, a detailed comparison of our final sample with the results of previous 
studies suggests that, at faint magnitudes, several high-redshift galaxy samples in the literature 
are significantly contaminated by low-redshift interlopers. 
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1 INTRODUCTION 

The goal of identifying and studying the nature of ultra high- 
redshift galaxies remains one of the most important challenges in 
observational cosmology, and holds the key to furthering our un- 
derstanding of the earliest stages of galaxy evolution and unveiling 
the nature of the sources responsible for cosmic reionization. 

Observational constraints provided by the Gunn-Peterson 
trough in the spectra of high-redshift quasars (e.g Fan et al. 2006) 
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suggest that reionization was coming to an end at z ~ 6 (Becker 
et al. 2007). Moreover, optical polarization measurements from the 
WMAP experiment indicate that reionization began at z ~ 1 1 if it 
is assumed to be a single, rapid event (Dunkley et al. 2009). Con- 
sequently, it is now apparent that to improve our understanding of 
cosmic reionization, and to unveil the earliest epoch of galaxy for- 
mation, it is necessary to extend studies of high-redshift galaxies 
into the redshift 7 < z < 10 regime (e.g. Robertson et al. 2010). 

Given our existing knowledge of the evolution of the galaxy 
luminosity function in the redshift interval 5.0 < z < 6.5 (e.g. 
Bouwens et al. 2007; McLure et al. 2009) it is clear that achieving 
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this aim requires ultra-deep near-IR imaging, reaching detection 
limits of 26 < tuab < 30. At the bright end of this range, wide- 
field, ground-based imaging has a unique contribution to make, and 
has recently allowed the luminosity function, the clustering prop- 
erties and the stellar populations of luminous Lyman-break galax- 
ies (LBGs: McLure et al. 2009; Ouchi et al. 2009; Grazian et al. 
2010) and Lyman alpha emitters (LAEs: Ouchi et al. 2010; Ota et 
al. 2010; Ono et al. 2010; Nakamura et al. 2010) to be studied in 
detail. Indeed, the importance of ground-based imaging and spec- 
troscopy has recently been highlighted by the spectroscopic con- 
firmation of two LBGs at z — 7.01 and z = 7.11 by Vanzella et 
al. (2010). However, in advance of 30-m class ground-based tele- 
scopes, it is clear that routinely identifying and studying sub-L* 
galaxies at z > 7 is only possible using space-based imaging. 

Consequently, the unparalleled near-IR sensitivity provided by 
the new WFC3 camera, installed on the Hubble Space Telescope 
(HST) in late 2009, has proven to be a crucial breakthrough in high- 
redshift galaxy studies. Indeed, despite only covering an area of 
~ 4.5 sq. arcmin, the unprecedented depth (rriAB — 29, 5a) of the 
first tranche of WFC3/IR imaging of the Hubble Ultra Deep Field 
(HUDF; GO- 1 1 563) led to a raft of early science papers investigat- 
ing the number densities, luminosity functions, stellar masses and 
stellar populations of 6.5 < z < 8.5 galaxies (e.g. Bouwens et al. 
2010ab; Oesch et al. 2010; McLure et al. 2010; Yan et al. 2010; 
Labbe et al. 2010; Finkelstein et al. 2010; Bunker et al. 2010). In- 
terestingly, Lehnert et al. (2010) have recently claimed the tenta- 
tive detection of Lya emission at z = 8.56 in a WFC3/IR candi- 
date in the HUDF, originally identified by Bouwens et al. (2010a) 
and McLure et al. (2010). Although the large rest-frame equiva- 
lent width (EW~ 210A) of the Lya emission line suggests that, 
if confirmed, this object must be a decidedly atypical example of a 
z > 6 LBG (Stark, Ellis & Ouchi 2010), the location of the claimed 
Lya emission line is in good agreement with the original photomet- 
ric redshifts derived by McLure et al. (2010) and Finkelstein et al. 
(2010); z phot = 8.45 ±0.50 and z phot = 8.61 ±0.35 respectively. 

In addition to the WFC3/IR imaging of the HUDF, another 
key dataset has been the WFC3/IR imaging taken as part of the 
Early Release Science extra-galactic programme (ERS; GO-1 1359) 
which, although substantially shallower than the HUDF WFC3/IR 
imaging (rriAB — 27.5, 5a), covers an area approximately ten 
times larger. By combining the HUDF and ERS datasets to obtain 
greater dynamic range in UV luminosity, Bouwens et al. (2010b) 
and Labbe et al. (2010) investigated the relationship between the 
UV spectral slope (fx oc X 13 ) and UV luminosity. Both studies 
find a correlation, with /3 changing from /3 ~ — 2 (a typical value 
for lower-redshift starburst galaxies) at M1500 — —20.5, to ex- 
tremely blue values of /3 ~ —3 at M1500 — —18.5. As dis- 
cussed by Bouwens et al. (2010b), although dust-free, low metal- 
licity models can produce slopes of j3 ~ —3, they can only do so 
under the assumption that the ionising photon escape fraction is 
high (f esc > 0.3) and that, correspondingly, the contribution from 
nebular continuum emission is low. 

Based on stacking the AC S +WFC3/IR+IR AC photometry of 
2 ~ 7 LBG candidates in the HUDF and ERS datasets, Labbe et al 
(2010) find the same correlation between Mi 500 and spectral slope 
/3 as Bouwens et al. (2010b). However, Labbe et al. (2010) conclude 
that it is not possible to reproduce both the blue spectral slopes 
and significant A rest ~ 4000A spectral breaks displayed by the 
faintest 2 ~ 7 LBG candidates, without recourse to episodic star- 
formation histories and/or a significant contribution from nebular 
line emission. Indeed, Ono et al. (2010) also conclude that nebular 



line emission may be necessary to reproduce the observed J— m 3 .6 
colour in a stack of 2 ~ 6 LAE photometry. 

Despite the uncertainties, one observational result that has re- 
ceived significant attention recently is the apparent relationship be- 
tween star-formation rate and stellar mass. Both Labbe et al. (2010) 
and Gonzalez et al. (2010) find an approximately linear correlation 
between stellar mass and star-formation rate at z ~ 7, consistent 
with the results derived by Stark et al. (2009) for LBGs in the red- 
shift range 4 < 2 < 6. As a result, Labbe et al (2010) and Gonzalez 
et al. (2010) conclude that the specific star-formation rate (sSFR) 
of 2 ~ 7 LBGs is remarkably constant (sSFR ~ 2 Gyr~ ), and 
consistent with the value of sSFR ~ 2.5 — 4.5 Gyr~ observed in 
star-forming galaxies at 2 ~ 2 — 3 (e.g. Daddi et al. 2007; Magdis et 
al. 2011). As previously discussed by Stark et al. (2009), a natural 
explanation of this observation would be to invoke a star-formation 
rate which exponentially increases with time although, as shown 
by Finlator, Oppenheimer & Dave (2011), star-formation histories 
of this type may have difficulty in reproducing some of the most 
extreme Balmer breaks reported in the literature at 2 > 6. 

The majority of previous studies which have investigated the 
high-redshift galaxy population using WFC3/IR imaging have re- 
lied on traditional colour-cut, or "drop-out", selection techniques. 
In contrast, the principal motivation for this paper is to investigate 
what can be learned about the 2 > 6 galaxy population by fully 
exploiting the excellent multi-wavelength (ACS+WFC3/IR+IRAC) 
data which is now available over an area of ~ 45 sq. arcmin. Rather 
than applying standard "drop-out" criteria, in this work we continue 
to pursue the strategy we have previously adopted (McLure et al. 
2006; 2009; 2010) and employ a template-fitting, photometric red- 
shift analysis to select our final high-redshift galaxy sample. A key 
new element in this strategy is our development of a deconfusion 
algorithm capable of providing the robust IRAC photometry neces- 
sary for improved photometric redshift and stellar-mass estimates. 

In principle this technique should have several advantages 
over the standard LBG "drop-out" selection. Firstly, by employ- 
ing all of the available multi-wavelength data, including the IRAC 
photometry, it is possible to make optimal use of the available in- 
formation. Secondly, by avoiding any colour pre-selection this ap- 
proach should be less biased towards simply selecting the very 
bluest galaxies at high-redshift. This second point is potentially 
crucial within the context of investigating the claims of ultra-blue 
UV spectral slopes for LBG candidates at 2 > 6 (see Section 4). Fi- 
nally, an SED fitting analysis also provides an estimate of the pho- 
tometric redshift probability density function, P(z), and therefore 
allows the prevalence, and significance, of competing photometric 
redshift solutions at low redshift to be transparently investigated. 

The primary motivation of this paper is therefore to construct 
the most robust sample possible using the available data and tech- 
niques reviewed above, in order to critically address some of the 
claims about the properties of the 2 > 6 population newly-found 
with HST. The structure of the paper is as follows. In Section 2 
we describe the available data in each of the three fields, including 
a brief description of our IRAC deconfusion algorithm. In Section 
3 we describe our initial candidate selection, photometric redshift 
analysis and the construction of the final catalogue of robust can- 
didates. In Sections 4 and 5 we investigate the UV spectral slopes, 
stellar masses and star-formation rates of the final robust sample. 
In Section 6 we perform a detailed comparison of our final robust 
sample with samples previously derived in the literature, explor- 
ing the reasons behind any apparent discrepancies. In Section 7 we 
provide a summary of our main conclusions. In the Appendix we 
provide a full description of our IRAC deconfusion procedure, full 
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photometry and grey-scale postage stamp images for each high- 
redshift candidate and individual plots illustrating the results of our 
SED fitting. All magnitudes are quoted in the AB system (Oke & 
Gunn 1983) and all calculations assume Qq — 0.3, Qa = 0.7 and 
H = rOkms^Mpc" 1 . 



2 DATA 

The analysis in this paper relies on the publicly available optical, 
near-IR and mid-IR imaging data covering the HUDF, HUDF09-2 
and ERS fields. In this section we briefly describe the relevant de- 
tails of the various imaging datasets, and the depth analysis which 
was performed in order to attribute accurate error estimate s to the 
candidate photometry. The basic properties of the three fields are 
listed in Table 1 . 

2.1 WFC3/IR imaging 

The WFC3/IR imaging of both the HUDF and the HUDF09-2 fields 
was taken as part of the public treasury programme GO- 11563 
(P.I.=niingworthfl and consists of single pointings (~ 4.5 sq. ar- 
cmin) of the WFC3/IR camera in the F105W, F125W and F160W 
filters (hereafter referred to as Y105 , J125 and Hieo). The WFC3/IR 
dataset in the ERS field was taken as part of the public programme 
GO-11359 (P.I.=0'Connell) and consists of a mosaic of 10 point- 
ings of the WFC3/IR camera in the F098M (Y 09 $), J 125 and H 160 
filters (Windhorst et al. 2011) 0. The WFC3/IR data were cali- 
brated using CALWF3 and subsequently combined using MultiDriz- 
zle (Koekemoer et al. 2002) as summarized in McLure et al. 2010); 
full details are presented in Koekemoer et al. (201 1). The final mo- 
saics have PSFs with FWHM in the range 0.15"— 0.1 8" depending 
on the filter, and were drizzled onto a final grid of 0.06"/pix. In 
the case of the HUDF and ERS mosaics, the final astrometry was 
matched to that of the publicly available reductions of the optical 
ACS imaging of the UDF (Beckwith et al. 2006) and GOODS-S 
(GOODSv2.0; Giavalisco et al. 2004) respectively. The astrome- 
try for the final mosaics of the HUDF09-2 field was matched to the 
J— band imaging of GOODS-S taken as part of the MUS YC survey 
(Cardamone et al. 2010), with a typical r.m.s. accuracy of ~ 0.1". 
The WFC3/IR imaging of the ERS analysed in this paper consists 
of the data comprising the completed programme. However, for the 
HUDF09-2 and HUDF fields we make use of the epoch 1 obser- 
vations, which consist of the data publicly available as of February 
2010 and August 2010 respectively. 

2.2 ACS imaging 

For the HUDF and ERS fields the ACS data used in this study con- 
sists of the publicly available reductions of the F435W, F606W, 
F775W and F850LP (hereafter -B435, Vf;o6, *775 and zsso) imag- 
ing of the HUDF (Beckwith et al. 2006) and GOODS-S (v2.0; Gi- 
avalisco et al. 2004). The ACS imaging covering the HUDF09-2 
field is our own reduction (based on CALACS and MultiDrizzle) 
of the Vsos, 2775 and zsso imaging obtained as part of the UDF05 

1 We do not consider the third WFC3/IR pointing obtained as part of GO- 
1 1563, HUDF09-1, because deep IRAC imaging of this field is not currently 
available. 

2 All of the WFC3/IR data utilised in this paper conforms to the nomi- 
nal flight zeropoints, i.e. Y098M = 25.68, Y W5W = 26.27, J125W = 
26.25, H 160W = 25.96. 



programme (Oesch et al. 2007). All of the optical ACS imaging 
was re-sampled to a 0.06"/pix grid to match the WFC3/IR data, 
and the astrometry of the ACS data covering HUDF09-2 was also 
registered to match the MUSYC imaging of GOODS-S. 



2.3 IRAC imaging 

For the HUDF and ERS fields we make use of the publicly avail- 
able reductions (v0.30) of the 3.6/itn and 4.5/im IRAC imaging ob- 
tained as part of the GOODS survey (proposal ID 194, Dickinson 
et al.; in preparation). The IRAC data covering the ERS consists of 
approximately 23 hours of on-source integration in both the 3.6/im 
and 4.5/im bands. For the HUDF we performed an inverse vari- 
ance weighted stack of the overlapping region of the epoch 1 and 
epoch 2 imaging, producing final mosaics consisting of approx- 
imately 46 hours of on-source integration at 3.6/im and 4.5/im. 
For the HUDF09-2 we re-registered and stacked the MOPEX reduc- 
tions of the 3.6/im imaging obtained via proposal ID 30866 (P.I. 
R. Bouwens) producing a final 3.6/im mosaic with an on-source 
integration time of approximately 33 hours. The consistency of 
the IRAC photometry was checked via reference to the SIMPLE 
(Damen et al. 2010) imaging at 3.6/im and 4.5/im of the E-CDFS, 
which overlaps all three fields. The astrometry of the IRAC imag- 
ing in all three fields was registered to that of the corresponding 
.Hieo WFC3/IR imaging. 



2.3.1 IRAC deconfusion 

A key new feature of the analysis undertaken in this paper is the 
inclusion of the IRAC photometry in the candidate high-redshift 
galaxy selection procedure. Due to its depth, and comparatively 
broad PSF (FWHM ~ 1.5"), the 3.6 + 4.5yum IRAC imaging 
covering the three fields of interest is heavily confused, making 
the process of obtaining aperture photometry matched to the opti- 
cal/nearlR HST imaging non trivial. In order to achieve this aim it 
is therefore necessary to pursue some form of deconfusion process, 
which allows the IRAC imaging to be utilised beyond the natural 
confusion limit. Although there are several techniques which can 
be used to deconfuse IRAC imaging (see Appendix A) we have de- 
veloped our own software which uses the WFC3/IR imaging data to 
provide normalized templates for each object in the field and then, 
via a transfer function, produces synthetic IRAC images on the na- 
tive 0.6"/pix plate scale. Through a matrix inversion procedure, the 
amplitude (or total flux) of each template can be simultaneously 
fitted to produce the optimal reproduction of the observed IRAC 
image (see Fig. 1). As a result of this procedure it is effectively 
possible to extract accurate aperture photometry from the 3.6/im 
and 4.5/im IRAC imaging at the spatial resolution of the WFC3/IR 
imaging. An additional advantage of this approach is that it natu- 
rally provides robust uncertainties on the delivered flux measure- 
ments, which depend both on the signal-to-noise ratio of the IRAC 
detection and the local level of confusion in the IRAC image. 



3 CANDIDATE SELECTION 

The process of candidate selection can be broken down into three 
separate stages: object detection and photometry, photometric red- 
shift analysis and sample cleaning. Each stage in the process is de- 
scribed below. 
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Figure 1. Illustration of the IRAC deconfusion algorithm. The left-hand panel shows the inverse-variance weighted stack of the epoch l+epoch2 4.5/im 
imaging covering the HUDF. The middle panel shows the best-fitting model of the IRAC data, based on using the -f/ieo WFC3/IR imaging to provide model 
templates, and a matrix inversion procedure to determine the best-fitting template amplitudes (see Appendix A for full details). The right-hand panel shows 
the model subtracted image (note that the WFC3/IR imaging does not cover the full area of the HUDF). 



Table 1. The results of the image depth analysis described in Section 3.2. Columns 1-4 list the survey field names, central coordinates and areas (in sq. 
arcmin). The areas listed in column 4 correspond to the areas of WFC3/IR imaging which were actually useful for the analysis performed in this paper 
(i.e excluding array edges and areas which are not covered by ACS optical imaging). The remaining columns list the median 5cr depths in each available 
ACS+WFC3/IR+IRAC filter. The depths are all referenced to a 0.6" diameter aperture, and are not aperture corrected. Due to the fact that the IRAC depths are 
determined via a deconfusion process (see Section 2.3.1), the local IRAC depths can differ significantly from the median values listed here. 



Field 


RA(J2000) 


DEC(J2000) 


Area 


B435 


^606 


*775 


^850 


^098 


Y W5 


J125 


#160 


3.6/im 


4.5/im 


HUDF 


03:32:38.5 


-27:46:57.0 


4.5 


29.04 


29.52 


29.19 


28.54 




28.59 


28.67 


28.73 


26.3 


25.9 


HUDF09-2 


03:32:23.4 


-27:42:52.0 


4.5 




28.49 


28.22 


28.06 




28.24 


28.60 


28.49 


26.2 




ERS 


03:33:05.5 


-27:51:21.6 


36.5 


27.68 


27.87 


27.29 


27.06 


27.26 




27.66 


27.40 


26.0 


25.6 



3.1 Object detection and photometry 

The initial catalog construction process was identical in each of the 
three fields, and relied on SEXTRACTOR v2.5.0 (Bertin et al. 1996). 
Preliminary catalogues were constructed in which object detection 
was performed in the Y\oz,w /Y09&M, J125VK and H\qqw bands, 
using an aggressive set of SEXTRACTOR parameters, with matched 
photometry extracted from the corresponding ACS imaging by run- 
ning SEXTRACTOR in dual image mode. The separate catalogues 
were then concatenated to produce a master catalog of unique ob- 
jects in each of the three fields. 

In order to avoid biases which can be introduced by adopt- 
ing small photometric apertures, all of the analysis in this paper 
is based on 0.6"diameter aperture photometry. For the purposes of 
the photometric redshift analysis, the fluxes from the 0.6"diameter 
aperture photometry are not corrected to total, but the WFC3/IR 
and IRAC fluxes are corrected by small amounts (2%-10%) to ac- 
count for aperture losses relative to the ACS imaging. 

3.2 Depth Analysis 

A crucial part of the analysis necessary to identify robust high- 
redshift candidates is the derivation of accurate photometric un- 
certainties in each band. This was achieved by first producing a 
so-called \ 2 image (Szalay, Connolly & Szokoly 1999) of the reg- 
istered optical+nearlR images of each field to identify which pixels 
are genuine "blank sky". Secondly, a grid of 0.6"diameter apertures 
was placed in the blank sky regions on each image. Thirdly, in or- 
der to determine the local image depth for each candidate, in each 
filter, the r.m.s. aperture-to-aperture variation was determined, by 
examining the distribution of the nearest 50 blank apertures. In this 
fashion, we are able to determine a local depth measurement for 
each individual candidate. For information, the median 5a depths 
for each field are listed in Table 1. 



3.3 Photometric redshift analysis 



To perform the SED analysis necessary for this study we have de- 
veloped a new, bespoke, template-fitting code. The primary moti- 
vation for developing this new code was to provide the freedom to 
explore the relevant multi-dimensional parameter space in detail, 
investigating the impact of different SED templates, IMFs, dust at- 
tenuation prescriptions and IGM absorption recipes. Moreover, by 
employing our own software it is possible to have full control over 
which derived quantities are provided as output, and the exact de- 
tails of how the template fitting is performed. For example, our new 
code performs the SED fitting based on flux densities (/„), rather 
than magnitudes, which has the advantage of allowing the flux er- 
rors to be dealt with in a rigorous manner. Moreover, if necessary, 
the new code offers the possibility of fitting the input photome- 
try with multi-component stellar populations, each with separate 
metallicities and/or dust attenuation prescriptions. 

For the purposes of this study we employed the Bruzual 
& Chariot (2003) and Chariot & Bruzual 2007 (priv. coram) 
stellar evolution models (hereafter BC03 & CB07), consid- 
ering models with metallicities ranging from solar (Zq) to 
l/50th solar (0.02Zq). Models with instantaneous bursts of star- 
formation, constant star-formation and star-formation rates ex- 
ponentially declining with characteristic timescales in the range 
50 Myrs < r < 10 Gyrs were all considered. The ages of the 
stellar population models were allowed to range from 10 Myrs to 
13.7 Gyrs, but were required to be less than the age of the Universe 
at each redshift. Dust reddening was described by the Calzetti et 
al. (2003) attenuation law, and allowed to vary within the range 
0.0 < Av < 2.5 magnitudes. Inter-galactic medium absorption 
short-ward of Lye* was described by the Madau et al. (1995) pre- 
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Figure 2. Example plots showing the results of our SED-fitting for three 
objects from our final robust sample at redshifts 2 p ^ ot = 6.5, 7.0 & 8.7 
(top to bottom). In each plot the thick (blue) line shows the best-fitting high- 
redshift galaxy template, and the thin dotted (red) line shows the best-fitting 
alternative solution at low redshift. In each panel the inset shows the value 
of x 2 as a function of redshift (marginalized over all other free parameters). 
The horizontal line in each inset panel highlights the location of (Xmi n + 4) 
which corresponds to our requirement that a candidate is considered "ro- 
bust" only if any alternative low-redshift solutions can be ruled-out at the 
> 95% confidence level. In each example, the combination of image depth 
and wavelength coverage allows us to rale-out any alternative low-redshift 
solutions with high confidence. The bottom panel (HUDF09-2_2765) pro- 
vides a good example of a situation where even a robust upper limit to the 
IRAC 3.6/^m flux provides a powerful redshift discriminant. 



scription, and a Chabrier (2003) IMF was assumed in all cases Bin 
Fig. 2 we show example SED fits for three objects (one from each 
field) covering the redshift range 6.5 < z p hot < 8.7. 



3.4 Sample cleaning 

Based on the results of the photometric redshift fitting, all objects 
which displayed a statistically acceptable solution at z p hot > 4.5 



3 Derived quantities such as stellar masses and star-formation rates can be 
converted to a Salpeter (1955) IMF by multiplying by a factor of 1.8. 
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Figure 3. In black we show the redshift distribution of our final robust sam- 
ple, which has been calculated by summing the estimated redshift proba- 
bility density function of each candidate. This is our best estimate of the 
redshift distribution of the N=70 objects which satisfy all three of the crite- 
ria listed in Section 3.4.1. In grey we show the redshift distribution of the 
N=130 objects which satisfy the first criterion listed in Section 3.4.1 (i.e. a 
statistically acceptable redshift solution at 2 p ^ot > 6.0) but fail the other 
two criteria (see text for details). 



were retained, while those with no acceptable solution at high red- 
shift were excluded. In each of the three fields, this initial screen- 
ing process removed more than 90% of the original input cata- 
logues. Following the initial photometric redshift fitting, the re- 
maining samples of potential high-redshift candidates were man- 
ually screened to remove artefacts (e.g. diffraction spikes), edge 
effects and spurious candidates such as high-surface brightness 
features within the extended envelopes of luminous low-redshift 
galaxies. 



3.4.1 Final candidate sample 

From a practical perspective, the primary goal of this study is 
to produce a robust sample of high-redshift galaxy candidates at 
z > 6. In order to achieve this aim, three criteria were applied to 
the remaining potential high-redshift candidates: 

• Statistically acceptable redshift solution at z p hot > 6.0 

• Secondary redshift solution excluded at > 95% confidence 

• Integrated probability J*Zq° P(z')Sz' > 0.5 

The first criterion simply restricts our final sample to those objects 
for which the best-fitting SED template lies at z p hot > 6. The sec- 
ond criterion rejects those objects for which the competing low- 
redshift solution cannot be excluded at high confidence. Specifi- 
cally, this criterion is enforced by insisting that the A\ 2 between 
the primary and secondary photometric redshift solution (following 
marginalization over all other relevant parameters) is A\ 2 > 4. 
The final criterion is designed to exclude a small number of can- 
didates with relatively flat P(z) distributions for which, despite 
having a primary photometric redshift solution at z p hot > 6, the 
majority of their integrated probability density distribution function 
lies at z p hot < 6. We note that this final criterion is very similar to 
that employed by Finkelstein et al. (2010) in their HUDF analysis. 
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Our final robust sample of z > 6 galaxy candidates consists 
of N=70 objects, spanning the redshift range 6.0 < z < 8.7 and 
covering more than a factor of ten in intrinsic UV luminosity, from 
— 18.2 < M1500 < —21.20. It is perhaps worth pointing out that 
if we had only insisted on a statistically acceptable primary pho- 
tometric redshift solution at z p hot > 6.0, the final sample would 
contain N=130 candidates. It should be stressed that it is likely that 
a significant fraction of the excluded objects are indeed z > 6 
galaxies (see Fig. 3), it is simply that with the data in-hand, it is not 
possible to consider them as robust candidates. 

Due to the fact that LBGs at z > 6 are necessarily young 
galaxies, the differences between the SED fits provided by the 
BC03 and CB07 models are negligible. Consequently, to ease the 
comparison with previous studies, throughout the rest of the pa- 
per we adopt the results of the SED-fitting analysis based on the 
BC03 models. The final robust samples in each of the three fields, 
along with the best-fitting photometric redshift solutions and vari- 
ous other derived parameters are presented in Tables 2, 3 & 4. Note 
that, although we make no further use of this information through- 
out the rest of the paper, in Tables 2, 3 & 4 we also list the best- 
fitting photometric redshift for each candidate if Lya emission with 
a rest-frame EWo in the range OA < EWo < 240A is included as 
an extra free parameter in the SED-fitting procedure. This infor- 
mation is provided to indicate the maximum plausible redshift for 
each candidate. The three candidates which are listed separately at 
the bottom of Tables 3 & 4 only pass our criteria as robust z > 6 
candidates if Lya emission is included in the SED-fitting proce- 
dure, and are not included in any of the subsequent analysis. The 
0.6"diameter aperture photometry for each candidate is listed in Ta- 
bles B 1 , B2 & B3 of Appendix B, along with plots of the best-fitting 
SED templates and grey-scale optical-nearlR postage stamps. 



4 THE UV SPECTRAL SLOPES 

As discussed in the introduction, one of the most interesting, and 
controversial, results to emerge from the new WFC3/IR-selected 
LBGs has been the claim that faint LBGs (A/1500 — — 18.5) at 
2 > 6 display extremely blue {J3 ~ —3) UV spectral slopes (e.g. 
Bouwens et al. 2010b; Labbe et al. 2010). Given the relatively small 
areas which have currently been imaged with WFC3/IR (i.e. ~ 50 
sq. arcmin), the brightest WFC3/IR-selected 2 ~ 7 LBGs have 
absolute UV luminosities of A/1500 — —21. It is widely agreed 
in the literature that at these absolute magnitudes (L ~ 2L*), 
2 ~ 7 LBGs display the same UV spectral slopes (/? ~ —2) as 
observed for young (~ 100 Myr) starbursts at redshifts 3 < z < 5. 
However, in contrast, it has been claimed that the faintest LBGs 
at 2 ~ 7 (Misoo — —18.5) display much bluer spectral slopes; 
(P) = -3.0 ± 0.2 (Bouwens et al. 2010b). 

Although this may appear to be a relatively small difference 
in spectral slope, it is potentially of great interest. The reason is 
very straightforward. While UV spectral slopes of f} ~ —2 can be 
comfortably reproduced by standard simple stellar population mod- 
els (without recourse to ultra-young ages or ultra-low metallicities), 
spectral slopes of /3 ~ — 3 cannot, and probably require a combina- 
tion of zero reddening, very young ages (i.e. ~ 10 — 30 Myrs) and 
a high escape fraction of photons short- ward of Lya (e.g. Bouwens 
et al. 2010b; Labbe et al. 2010). Given the potential importance of 

4 The final sample consists of 73 objects if 3 additional objects are included 
which only satisfy our selection criteria when Lya emission is included in 
the SED templates. 
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Figure 4. Plots of UV spectral slope (J3) versus redshift for the final robust 
sample of seventy objects at z p hot ^ 6. The top three panels show ft versus 
redshift for the three separate survey fields, while the bottom panel shows 
ft versus redshift for the full combined sample. In each panel the horizontal 
line shows the variance-weighted mean value of ft (see Table 5). 



this result, not least for studies of reionization, it is clearly of inter- 
est to investigate the UV spectral slopes displayed by the sample of 
high-redshift LBGs derived here. 

The individual values of /3 measured for each candidate are 
listed in Tables 2, 3 & 4. The f3 values have been calculated using 
the following formulae: 
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Table 2. Details of the final high-redshift sample in the HUDF field. The first three columns list the ID number and coordinates of each candidate. Columns 
four, five and six list the photometric redshifts, the lcr uncertainty on the photometric redshifts and the values of x f° r the best-fitting SED template. Column 
seven lists the difference in \ 2 between the best-fitting high-redshift SED template and the alternative photometric redshift solution at low redshift (typically 
at 1.0 < Zphoi < 2.5). Column eight lists the absolute magnitude of the best-fitting SED template, where M1500 is measured using a lOOA-wide filter 
centred on a rest-frame wavelength of 150()A. Column nine lists the value of the UV spectral slope (fx oc A' 3 ) for each candidate, derived using the formulae 
listed in Section 4. Column ten lists the best-fitting photometric redshift derived including Lya emission as an additional free parameter, and column eleven 
lists the logarithm of luminosity of the corresponding best-fitting Lya emission line (in ergs s~ x ). Column twelve lists the results of cross-checking each 
candidate against existing literature studies of the WFC3/IR data in the HUDF field. Matches were found with the following papers: M (McLure et al. 2010); 
B (Bouwens et al. 201 1; robust), Bp (Bouwens et al. 2011; potential), F (Finkelstein et al. 2010), W (Wilkins et al. 2010), Y (Yan et al. 2010) and L (Lorenzoni 
etal. 2011). 



ID 


RA(J2000) 


DEC(J2000) 


z 


A 


2 


x 2 


A X 2 


M1500 


P 




\og{L Lya ) 


Literature 


UT TDF 1 ^44 


UJ. JZ.JD.DJ 


_27.47.50 [ 


6 06 


5 89- 


-6 17 


1 1 


27 6 


1 Q ^ 


9 9 4-0 7 
— z.z u. 1 


6 16 


42 2 


M 


HTTDF 1016 


03-39-35 ()f\ 

U J .Ji 1 J ,IAJ 


— 27'47'40 2 


6 06 


5 97- 


-6 15 


4 3 


32 


— 19.6 


—0.4 ± 0.6 


6 33 


42.4 


M 


HTTDF S99 


n3-39-3fi ai 


97-Af,-A\ zl 
— z 1 .40.4 1 .4 


A 07 
o.u / 


J.70- 


U. 1 J 


zl 5 


iji.j 


90 R 


1 Q4-(lfi 








HTTDF IfOl 




97-47-SO 9 

— Z 1 .4/ JU.Z 


f\ 1 1 


J.yJ 


A 38 

O.JO 


1 ^ 
1 .j 


1 ^ zt 

1J4 


1 8 R 


2 2 ± 1 2 


6 4^ 


49 ^ 


M 


HTTDF 7Q6 
nuur./w 




97-zt6-19 8 

— Z ( .4U.JZ.O 


f\ 1 Q 
0.17 


S SA 
J.oU 


A 31 
O.J 1 


1 4 


a& n 

40.VJ 


in n 

— ly.y 


1 7 4- 
1 . / x u.u 


6 SO 


49 7 


M F 


HTTDF 78^6 


n3-39-35 05 


97-Z17-9S 8 

— Z( .4/ .ZJ.O 


6 22 


5 Q7 
J.y 1 


U.4J 


a 

U.O 


7 A 


^g 4 


A f\ 4- 9 
— 4.u in z.o 


6 SI 

O.J 1 


42 1 


M 


HTTDF 16Q9 


03-39-43 D3 


— 97-46-93 6 

Z 1 .4U.ZJ .U 


6.23 


6.1 1- 


-6.34 


3.1 


38.5 


— 19.3 


—3.4 it 0.8 


6.49 


42.4 


M 


HTTDF 11 'Al 


n3'39-3fi 


— Z J .4U.4Z.U 


6 26 


5 80- 


-6 72 


U.J 


4 


lg 2 


— 1 2 it 2 


6 40 


41 7 


M Y 


HTTDF 9^1 6 


UJ.JZ.44.J 1 


97-ztA-zlS 9 

— Z J .4D.4J.Z 


6 31 


6 03- 


-6 54 


1 2 


Q 1 

y.i 


18 7 


17 4-17 


6 ^0 




M 


HTTDF 99R1 


Uj.jZ.j7. / y 


_ 27 - 46'33 7 


6 37 


6 11- 


-6 57 


3 


5 7 


lg g 


q 9 4_ 1 7 
— o.z m j.. 1 


6 35 




M 


HTTDF 144? 


uj.jz.4i. ly 


97-ZIA-97 8 

— Z J .4U.Z / . O 


f\ 37 

U.J / 


6 17- 


-6 55 


6 4 


112 


10 2 


1 A 4- O 7 


6 43 


41 7 


ivi,r% vv,_ 


HTTDF 7194 


03-32-41 AO 

UJ . Ji.l 1 .KJ\J 


—27:47:04.5 


6.41 


6.18- 


6.60 


0.6 


5.5 


— 18.6 


—2.9 it 1.4 


6.40 




B 


HUDF.2672 


03:32:37.80 


-27:47:40.4 


6.45 


6.14- 


-6.67 


0.5 


7.9 


-18.6 


-2.0 it 1.2 


6.81 


42.4 


M 


HUDF.1818 


03:32:36.38 


-27:47:16.3 


6.57 


6.35- 


-6.72 


2.1 


17.3 


-19.1 


-2.3 it 1.0 


7.05 


42.7 


M,F,W,B,Y 


HUDF.1473 


03:32:36.77 


-27:47:53.6 


6.57 


6.42- 


-6.71 


2.1 


24.8 


-19.2 


-1.8 it 0.5 


6.99 


42.6 


M,F,W,B 


HUDF.1730 


03:32:43.78 


-27:46:33.7 


6.60 


6.37- 


6.84 


0.5 


11.5 


-19.1 


-2.7 it 0.8 


6.59 




M,F,W 


HUDF.1632 


03:32:37.44 


-27:46:51.2 


6.60 


6.40- 


-6.74 


0.7 


13.3 


-19.1 


-2.1 ±0.8 


6.60 




M,F,W,B,Y 


HUDF.2084 


03:32:40.57 


-27:46:43.6 


6.61 


6.39- 


-6.80 


2.6 


11.9 


-18.8 


-2.2 ±0.9 


7.13 


42.6 


M,F,W,B,Y 


HUDF.1995 


03:32:39.58 


-27:46:56.5 


6.62 


6.31- 


-6.91 


4.2 


6.2 


-18.9 


-3.5 ±0.9 


6.60 




M,F,B,Y 


HUDF.658 


03:32:42.56 


-27:46:56.6 


6.63 


6.53- 


6.79 


1.4 


81.9 


-20.5 


-2.0 ±0.5 


6.85 


42.5 


M,F,W,B,Y 


HUDF.2701 


03:32:41.82 


-27:46:11.3 


6.66 


6.35- 


-6.91 


2.3 


4.5 


-18.5 


-3.3 ±1.5 


6.88 


42.0 


F,W,Bp,Y 


HUDF.860 


03:32:38.81 


-27:47:07.2 


6.96 


6.72- 


-7.23 


1.8 


31.8 


-20.0 


-1.4 ±0.5 


6.96 




M,F,W,B,Y 


HUDF.1102 


03:32:39.55 


-27:47:17.5 


7.06 


6.75- 


-7.42 


2.5 


7.1 


-19.7 


-1.5 ±0.5 


7.06 




M,F,B,Y 


HUDF.1419 


03:32:43.13 


-27:46:28.5 


7.23 


6.80- 


-7.48 


5.9 


9.1 


-19.2 


-1.3±0.7 


7.95 


42.8 


M,F,W,B,L,Y 


HUDF.2641 


03:32:39.73 


-27:46:21.3 


7.35 


6.97- 


-7.76 


1.2 


7.2 


-18.7 


-4.3 ±1.4 


8.06 


42.6 


M,F,B,Y 


HUDF.1962 


03:32:38.36 


-27:46:11.9 


7.36 


6.80- 


-7.73 


1.0 


5.5 


-19.1 


-3.2 ±1.1 


7.27 




B,F,Y 


HUDF.1173 


03:32:44.70 


-27:46:44.3 


7.36 


7.07- 


-7.72 


4.5 


9.0 


-19.8 


-2.4 ±0.6 


7.36 




M,F,B,Y 


HUDF.2664 


03:32:33.13 


-27:46:54.5 


7.45 


6.98- 


7.89 


1.9 


4.2 


-18.6 


-4.8 ±2.4 


8.08 


42.5 


M,B,L 


HUDF.1660 


03:32:37.21 


-27:48:06.2 


7.52 


7.24- 


-7.76 


0.9 


14.4 


-19.3 


-2.6 ±0.7 


7.98 


42.5 


M,F,B,Y 


HUDF.1679 


03:32:42.88 


-27:46:34.5 


7.88 


7.51- 


-8.11 


1.7 


5.5 


-19.1 


-2.0 ±0.8 


8.80 


42.7 


M,F,B,L,Y 


HUDF.2003 


03:32:38.13 


-27:45:54.0 


8.49 


8.08- 


-8.75 


0.9 


7.5 


-19.1 


-2.3 ±1.1 


8.89 


42.6 


M,F,B,L,Y 



p = 4.43(Ji2s - H 160 ) - 2.0 (1) 
p = 5.47(yioB - J125) - 2.0 (2) 

,3 = 3.91(^098 - J125) -2.0 (3) 

depending on the available filters and the redshift of the candidate. 
To derive the above formulae we have adopted the following pivot 
wavelengths for the 11)98, Vios, </i25 & #160 filters: 0.9864/im, 
1.0552^m, 1.2486/jim & 1.5369>m (WFC3 Instrument Handbook 
for Cycle 19). In order to sample as similar a range of rest-frame 
wavelengths as possible, and to ensure no potential contamination 
from Lya line emission, the values of P have been calculated using 
equations 2 or 3 for candidates at z p hot _ 6.4 and equation 1 for 
those candidates at z v hot > 6.4. In Fig. 4 we plot the estimated 
UV spectral slopes versus photometric redshift for the final robust 
sample, split by field. Several features of this plot are worthy of 
comment and are briefly discussed below. 



4.1 Uncertainties on derived UV spectral slopes 

As can readily be seen from Fig. 4, the uncertainties on measuring 
P are typically large. This is simply a consequence of attempting to 
determine a spectral slope using two filters which are not well sep- 
arated in wavelength. As an illustration, consider a galaxy at 2 = 7 
with a canonical UV spectral slope of P = —2, which is detected at 
5ct significance in both the J125 and ij?i6u filters. The correspond- 
ing estimate of the UV spectral slope is P = —2.0 ± 1.3, where the 
error simply reflects the photometric uncertainty. Clearly, deriving 
meaningful estimates of P on an individual object-by-object basis 
requires significantly better than 5a photometry in both filters. One 
obvious method of overcoming this problem is to assume that each 
P measurement, although inaccurate, is at least unbiased. In which 
case, one can proceed to bin the data and attempt to estimate the 
mean value of p. However, even when adopting this approach, it 
is necessary to account for the wide range in P uncertainties dis- 
played by the objects in a typical sample, by calculating a properly 
weighted mean: 
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Table 3. Details of the final high-redshift sample in the ERS field. Columns one to eleven list the same quantities as in Table 2. Column twelve lists the results 
of cross-checking each candidate against existing literature studies of the WFC3/IR data in the ERS field. Matches were found with the following papers: B 
(Bouwens et al. 2011; robust), Bp (Bouwens et al. 2011; potential), W (Wilkins et al. 2010) and L (Lorenzoni et al. 2011). 



ID 


RA(J2000) 


DEC(J2000) 


z 


Az 


2 

X 


A 2 


M1500 


P 


^Lya 


\og(L Lya ) 


Literature 


ERS.7086 


03:32:34.75 


-27:40:35.1 


6.18 


6.02-6.35 


1.4 


11.2 


-20.2 


-2.4 ± 0.6 


6.36 


42.5 




ERS.6066 


03:32:07.86 


-27:42:17.8 


6.19 


5.87-6.41 


2.7 


15.8 


-20.3 


-2.5 ± 0.6 


6.66 


43.1 




ERS.9869 


03:32:15.40 


-27:43:28.6 


6.21 


6.01-6.41 


0.4 


8.0 


-19.8 


-3.4 ± 1.0 


6.44 


42.5 


Bp 


ERS.8668 


03:32:27.96 


-27:41:19.0 


6.22 


5.88-6.55 


2.0 


7.6 


-19.9 


-2.2 ± 0.8 


6.22 






ERS.9100 


03:32:20.24 


-27:43:34.3 


6.27 


5.95-6.50 


1.4 


6.6 


-19.8 


-2.3 ± 0.6 


6.45 


42.5 


Bp 


ERS.7225 


03:32:36.31 


-27:40:15.0 


6.30 


6.02-6.73 


4.3 


11.4 


-20.1 


-1.7 ± 0.7 


7.16 


43.1 




ERS.6438 


03:32:25.28 


-27:43:24.2 


6.33 


6.14-6.70 


6.7 


9.5 


-20.3 


-2.0 ± 0.7 


7.26 


43.2 


W 


ERS.6263 


03:32:06.83 


-27:44:22.2 


6.36 


6.14-6.59 


4.9 


10.0 


-20.3 


-2.1 ± 0.8 


6.40 


41.9 


B 


ERS.7776 


03:32:03.77 


-27:44:54.4 


6.46 


6.15-6.66 


3.9 


6.5 


-20.0 


-1.6 ± 0.8 


6.58 


42.3 




ERS.5847 


03:32:16.00 


-27:43:01.4 


6.49 


6.31-6.59 


2.7 


17.2 


-20.5 


-2.4 ± 0.8 


6.90 


43.1 


W 


ERS.8987 


03:32:16.01 


-27:41:59.0 


6.52 


6.06-6.84 


2.5 


5.9 


-19.7 


-1.8 ± 1.0 


6.72 


42.5 


B 


ERS.3679 


03:32:22.66 


-27:43:00.7 


6.55 


6.42-6.71 


5.5 


14.5 


-21.2 


-1.9 ± 0.5 


6.55 




W,Bp 


ERS.7412 


03:32:09.85 


-27:43:24.0 


6.57 


6.37-6.77 


9.1 


8.8 


-20.2 


-0.7 ±0.7 


7.58 


43.3 




ERS.6427 


03:32:24.09 


-27:42:13.9 


6.65 


6.37-6.88 


1.2 


10.3 


-20.3 


-2.5 ± 0.6 


6.64 


42.2 


W,B 


ERS.8858 


03:32:16.19 


-27:41:49.8 


6.79 


6.33-7.08 


1.2 


6.8 


-20.0 


-2.9 ± 0.8 


6.77 




B 


ERS.7376 


03:32:29.54 


-27:42:04.5 


6.79 


6.50-6.98 


0.6 


5.3 


-20.2 


-2.4 ±1.1 


7.27 


43.0 


W,B 


ERS.8176 


03:32:23.15 


-27:42:04.7 


6.81 


6.62-6.98 


4.4 


13.6 


-20.1 


-3.8 ± 1.1 


7.73 


43.3 


W 


ERS.7672 


03:32:10.03 


-27:45:24.6 


6.88 


6.64-7.05 


3.8 


7.7 


-20.3 


-2.4 ±1.0 


7.77 


43.3 




ERS.7475 


03:32:32.81 


-27:42:38.5 


7.11 


6.83-7.31 


5.2 


6.2 


-20.3 


-1.9 ± 0.7 


7.76 


43.2 




ERS.7236 


03:32:11.51 


-27:45:17.1 


7.18 


6.99-7.35 


5.0 


5.4 


-20.3 


-3.4 ±0.8 


7.74 


43.0 




ERS.9041 


03:32:23.37 


-27:43:26.5 


8.02 


7.61-8.20 


6.2 


9.1 


-20.0 


-5.0 ± 1.5 


8.16 


43.0 


L 


ERS.10288 


03:32:35.44 


-27:41:32.7 


8.28 


7.59-8.56 


2.1 


7.8 


-20.1 


-2.3 ± 0.8 


9.50 


43.2 


B 


ERS.8584 


03:32:02.99 


-27:43:51.9 


8.35 


7.63-8.74 


2.9 


4.7 


-20.4 


-1.6 ± 1.1 


9.37 


43.3 


B,L 


ERS.8496 
ERS.9923 


03:32:29.69 
03:32:10.06 


-27:40:49.9 
-27:45:22.6 


6.07 
6.59 


5.49-6.56 
6.37-6.78 


4.6 
6.8 


9.0 
7.7 


-19.7 
-20.0 


-1.8 ± 1.1 
-1.7 ± 1.3 


6.87 
7.59 


42.9 
43.3 





Table 4. Details of the final high-redshift sample in the HUDF09-2 field. Columns one to eleven list the same quantities as in Table 2. Column twelve lists the 
results of cross-checking each candidate against existing literature studies of the WFC3/IR data in the HUDF09-2 field. Matches were found with the following 
papers: B (Bouwens et al. 2011; robust), Bp (Bouwens et al. 2011; potential) and W (Wilkins et al. 2010). 



ID 


RA(J2000) 


DEC(J2000) 


z 


A 


2 


x 2 


Ax 2 


Ml 500 


/3 


ZLya 


\og(L Lya ) 


Literature 


HUDF09-2.2459 


03:33:06.30 


-27:50:20.2 


6.06 


5.92- 


6.17 


3.2 


20.8 


-19.2 


-1.4 ±1.0 


6.37 


42.5 




HUDF09-2.2613 


03:33:06.52 


-27:50:34.6 


6.08 


5.90- 


-6.22 


1.0 


13.2 


-19.2 


-3.0 ± 1.0 


6.44 


42.5 




HUDF09-2.2638 


03:33:06.65 


-27:50:30.2 


6.14 


5.76- 


6.46 


0.1 


4.5 


-19.0 


-2.1 ± 1.1 


6.44 


42.2 




HUDF09-2.1543 


03:33:01.18 


-27:51:22.3 


6.18 


6.05- 


6.26 


0.3 


22.3 


-20.4 


-2.3 ±0.6 


6.12 






HUDF09-2.605 


03:33:01.95 


-27:52:03.2 


6.30 


6.06- 


6.49 


0.1 


6.0 


-19.5 


-1.9 ±0.7 


6.30 






HUDF09-2.2587 


03:33:04.20 


-27:50:31.3 


6.30 


6.11- 


6.39 


3.3 


27.9 


-20.3 


-1.3 ±0.6 


6.28 


41.8 




HUDF09-2.1660 


03:33:01.10 


-27:51:16.0 


6.36 


6.15- 


6.47 


3.8 


7.9 


-20.3 


-1.1 ±0.6 


6.26 






HUDF09-2.1745 


03:33:01.19 


-27:51:13.3 


6.52 


6.22- 


-6.82 


0.3 


7.9 


-19.4 


-1.9 ±0.6 


6.98 


42.7 


W,B 


HUDF09-2.1620 


03:33:05.40 


-27:51:18.8 


6.61 


6.30- 


6.93 


1.7 


7.0 


-19.1 


-1.7 ±1.0 


7.39 


42.9 


W,B 


HUDF09-2.1721 


03:33:01.17 


-27:51:13.9 


6.73 


6.39- 


-7.05 


2.5 


4.3 


-19.8 


-2.1 ±0.5 


6.78 


42.0 




HUDF09-2.2455 


03:33:09.65 


-27:50:50.8 


6.82 


6.73- 


6.89 


1.9 


34.4 


-20.6 


-1.9 ±0.5 


7.12 


43.1 


W,Bp 


HUDF09-2.1584 


03:33:03.79 


-27:51:20.4 


7.17 


6.79- 


7.36 


0.7 


15.6 


-20.6 


-1.6 ±0.5 


8.03 


43.3 


W,B 


HUDF09-2.2814 


03:33:07.05 


-27:50:55.5 


7.30 


6.90- 


7.66 


0.5 


5.7 


-19.7 


-2.5 ±0.6 


7.26 




Bp 


HUDF09-2.1596 


03:33:03.76 


-27:51:19.7 


7.45 


7.06- 


-7.62 


6.0 


15.8 


-20.4 


-2.0 ±0.5 


7.95 


43.1 


B 


HUDF09-2.2000 


03:33:04.64 


-27:50:53.0 


7.68 


7.30- 


7.90 


1.7 


15.5 


-19.7 


-2.2 ±0.5 


8.01 


42.7 


B 


HUDF09-2.2765 


03:33:07.58 


-27:50:55.0 


8.70 


8.37- 


9.05 


0.9 


6.9 


-20.0 


-0.5 ± 0.6 


8.76 


41.8 


B 


HUDF09-2.799 


03:33:09.15 


-27:51:55.4 


6.88 


6.70- 


-7.00 


9.1 


15.5 


-19.5 


-1.6 ±0.6 


7.67 


43.1 


B,W 



Pi where f3i represents an individual /3 measurement for a single can- 

cr 2 didate, and of is the corresponding variance. 

(?) = ^ir—r (4) 
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Table 5. Estimates of the typical value of the UV spectral slope (/J) for our 
final robust sample. The first two columns list the names and sizes of the dif- 
ferent samples being considered. Columns three and four list the variance- 
weighted and arithmetic mean values of /3 respectively, together with their 
corresponding uncertainties. Note that candidate HUDF09-2.2765 has been 
excluded from these calculations because it provides a biased estimate of /3 
due to its high redshift (z phot = 8.7 ± 0.3). 



Sample 


N 


\Pvar / 


</3arith> 


HUDF 


31 


-1.99 ±0.14 


-2.40 ±0.18 


HUDF09-2 


15 


-1.90 ±0.15 


-1.91 ±0.13 


ERS 


23 


-2.32 ±0.18 


-2.52 ±0.19 


COMBINED 


69 


-2.05 ±0.09 


-2.33 ±0.11 



4.2 Average values of U V spectral slopes 

The variance-weighted values of (/3) for each sub-sample, and the 
full combined sample are listed in Table 5 where, for comparison, 
we also list the straight arithmetic means and standard errors. It 
can be seen from Table 5 that for the HUDF and ERS sub-samples 
(and for the full combined sample), the variance-weighted mean 
results in a significantly redder estimate of the typical value of the 
UV slope than the straight arithmetic mean. Interestingly, for the 
HUDF09-2 sub-sample, where the photometry is most robust (see 
discussion below), the difference between the two estimates is neg- 
ligible. 

The results listed in Table 5 indicate that the ERS sub-sample 
contains a higher percentage of objects with j3 < — 2 than the 
other two fields (based on the variance-weighted means), although 
the difference is not significant. However, it is worth noting that 
any suggestion that the ERS candidates display bluer UV spectral 
slopes cannot be due to a trend for increasingly blue UV spectral 
slopes with decreasing UV luminosity, given that the median abso- 
lute magnitude of the ERS sample is M1500 = —20.2 compared 
to M1500 = —19.1 for the HUDF. Overall, our results provide no 
evidence that the members of the 2 > 6.5 LBG population display 
values of /3 significantly different from those seen in comparably 
luminous LBGs in the redshift interval 3 < z < 5. 

4.3 Potential for bias 

It can be seen from Fig. 4 that the HUDF09-2 sub-sample seems 
to display a particularly tight distribution of UV slopes, whereas 
the HUDF and ERS sub-samples show considerably more scatter. 
At least part of the explanation for this is that the HUDF09-2 sub- 
sample has the most robust WFC3/IR photometry. The reason is 
that, although the WFC3/IR imaging of HUDF09-2 is deep (par- 
ticularly the J125 data), the supporting data at other wavelengths is 
not, in a relative sense, as good (e.g. no B435 data, relatively shal- 
low V606+*77B data, and no 4.5pim data). As a consequence, candi- 
dates in the HUDF09-2 field are required to be somewhat brighter 
in the near-IR in order to pass our robustness criteria (see photom- 
etry in Appendix B). 

Another noteworthy point is that the bluer mean UV slope 
in the ERS is probably connected to the relative depths of the 
WFC3/IR imaging in this field. Due to the fact that the J125 imag- 
ing in the ERS is significantly deeper than the accompanying Yogs 
and Hieo imaging, the ERS sub-sample is the closest of the three to 
being purely J125 — selected. It is clear that when estimating the UV 
spectral slope from the J125 — iiieo colour, selecting the sample 
largely on the apparent J125 magnitude must introduce the poten- 



tial for biasing the sample towards objects with blue values of /?. A 
proper investigation of the sources of bias, and the potential for con- 
straining the true underlying distribution of UV spectral slopes, re- 
quires detailed simulation work which, although beyond the scope 
of this paper, is investigated in detail by Dunlop et al. (2011). 



5 STELLAR MASSES AND STAR-FORMATION RATES 

A key advantage of employing a template-fitting SED analysis is 
that stellar mass and star-formation rate estimates can be directly 
derived from the best-fitting models. In this Section we use this 
information to investigate the relationship between stellar mass 
and UV luminosity, and try to determine the typical specific star- 
formation rate (sSFR) for L* LBGs at 2 > 6. Throughout this 
Section we will repeatedly refer to the results for a sub-sample 
of twenty-one objects which have the most reliable star-formation 
rate and stellar-mass estimates because they are detected at either 
3.6/im or 3.6 + 4.5/im. In Table 6 we list the best-fitting parame- 
ters returned by our SED-fitting analysis for these objects, based 
on the best-fitting templates drawn from the full range of star- 
formation histories, metallicities and reddening described in Sec- 
tion 3.3. Based on photometry alone it is very difficult to accurately 
constrain the star-formation history and metallicity of high-redshift 
galaxies. As a result, it is common in the literature to derive star- 
formation rate and stellar-mass estimates from a much more re- 
stricted set of SED templates (typically constant star-formation rate 
models). In order to investigate the effect of this approach, in Ta- 
ble 6 we also list the best-fitting parameters derived from an SED 
template with a constant star-formation rate, Z = 0.2Zq and zero 
reddening (CSF model). 

5.1 Stellar mass - UV luminosity relation 

In Fig. 5 we plot stellar mass versus absolute UV magnitude 
(M1500) for our final robust sample of seventy 2 > 6 LBGs. In the 
left-hand panel we plot stellar-mass estimates based on the best- 
fitting SED templates drawn from the full grid of star-formation 
histories, metallicities and reddening. In the right-hand panel we 
plot the stellar-mass estimates based on the CSF model alone. In 
both panels the small open circles indicate those objects which are 
formally undetected at 3.6/im. For these objects the only stellar 
mass constraints at A rcst > 4000A come from the upper limits 
at 3.6/4.5^im provided by the deconfusion analysis. In contrast, 
those objects which are detected at 3.6/im (> la) are plotted as 
large grey circles, and those objects detected at both 3.6/im and 
4.5/im are plotted as large black circles. 

Based on the data presented in left-hand panel of Fig. 5 we 
have used the FITXY routine (Press et al. 1992) to derive the follow- 
ing relationship between stellar mass and UV luminosity (L1500): 

log I vr 1 I = (2.14 ±0.56) log ( Ll5 "" 1 | - 37.05 ± 4.52(5) 
V M ©/ V WHz / 

which is shown as the thick black line in the left-hand panel 
of Fig. 5. It is interesting to compare our equation 5 with the 
M* — Ltjv relation derived by Gonzalez et al. (2011) based on 
~ 500 B— drop galaxies at 2 ~ 4. The M t — Ltjv relation derived 
by Gonzalez et al. (2011) is plotted as the thick grey line in both 
panels of Fig. 5, and has the form: M* oc L 15 oo°' 2 . It can be seen 
from the left-hand panel of Fig. 5 that both relations are clearly 
consistent, although the M* — Ltjv relation derived here for LBGs 
with a mean redshift of 2 = 6.8 ± 0.1 is somewhat steeper. 
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Figure 5. Stellar mass versus absolute UV magnitude (M1500) for the final robust sample of seventy z > 6 LBGs, where M1500 i s derived from the best-fitting 
SED models using a 100 A - wide filter centred on a rest-frame wavelength of 1500 A. In each panel, the small open circles are objects with non-detections (i.e. 
< 2a) at 3.6/im, the large grey circles are objects which are detected at 3.6/im and the large black circles are objects detected at both 3.6/im and 4.5/im. In the 
left-hand panel the stellar-mass measurements are based on our SED-fitting analysis using the full range of star-formation histories, metallicities and reddening 
described in Section 3.3. In the right-hand panel the stellar-mass measurements are based on a single set of SED templates with a constant star-formation rate, 
Z = 0.2Zq and zero reddening (CSF model, see text for details). In the left-hand panel the thick black line is our best-fitting M* — Luv relation (Equation 
5) and the dotted line indicates the lower limit enforced by insisting that each candidate has an age > 10 Myrs. Those objects which lie on the lower limit 
(small grey open circles) have very poorly constraint stellar masses and were not included in the derivation of the best-fitting — L V y relation. In both 
panels the thick grey line is the M* — Luv relation (corrected to Chabrier IMF) derived by Gonzalez et al. (201 1) using a large sample of LBGs at z ~ 4. 
In the right-hand panel the upper dotted line indicates the expected M t — Luv relation for a CSF model which has been forming stars for ~ 600 Myrs (i.e. 
since z ~ 20 for an object at the mean redshift of the final robust sample; z = 6.8 ± 0.1). 



In an earlier study, Stark et al. (2009) also explored the 
AL» — Luv relation based on 4 < 2 < 6 LBGs selected from 
the GOODS N+S fields. At z ~ 4 the data from Stark et al. (2009), 
based on a sample of ~ 700 B— drop candidates, is entirely consis- 
tent with the M* — Luv relation derived by Gonzalez et al. (201 1). 
At higher redshifts, both Stark et al. (2009) and Gonzalez et al. 
(2011) investigated the M* — Luv relation at z ~ 5 and z ~ 6, 
based on samples of V — drop and i— drop LBG candidates respec- 
tively. Interestingly, at z ~ 5 the results of both studies do appear 
to be consistent with a steepening of the M* — Luv relation. In 
fact, this effect was noted by Stark et al. (2009) but, based on the 
available data, both authors concluded that there was no strong ev- 
idence for redshift evolution. At z ~ 6 neither study had sufficient 
dynamic range in Luv to constrain the slope of the AL, — Luv 
relation. 

At a given UV luminosity, the range of stellar masses dis- 
played by the LBG candidates in Fig. 5 is simply a function of their 
M*/Luv ratios which, in turn, are largely a function of their stellar 
population ages. Unfortunately, those candidates which have nei- 
ther detections or meaningful upper limits at IRAC wavelengths in- 
evitably have stellar ages/masses which are very poorly constrained 
(small grey open circles in Fig. 5). These objects (which we have 
excluded from our determination of the best-fitting AL* — Luv re- 
lation) can be seen to congregate close to the lower limit which 
is imposed during the SED-fitting procedure by insisting that each 
candidate must have an age of > 10 Myr. However, in reality, the 
majority of these objects can tolerate SED fits with stellar popula- 
tions as old as ~ 200 Myrs, at which point their estimated stellar 



masses become an order of magnitude larger. Consequently, the 
apparent steepening of the M* — Luv relation at faint magnitudes 
must be viewed with considerable caution. 

It is clear from Fig. 5 that, based on the current sample, it 
is not possible to determine if the AL, — Luv relation at z > 6 
is steeper than at z ~ 4. Indeed, our results for the twenty- 
one objects with the most reliable stellar-mass estimates are en- 
tirely consistent with the conclusion that the slope and normali- 
sation of the AL — Luv relation does not change over the red- 
shift interval 4.0 < z < 7.0. However, by restricting ourselves 
to those objects with the most reliable stellar-mass estimates, the 
results presented in the left-hand panel of Fig. 5 suggest that L* 
(M1500 — —20.2) galaxies at z ~ 6.5 have a median stellar mass 
ofAL, = (2.1 ± 1.1) x 10 9 M . Moreover, by deriving stellar- 
mass estimates using stellar population models covering a wide 
range of metallicities, star-formation histories and reddening, our 
results indicate that the full range of M*/Luv ratios displayed by 
L* galaxies at this epoch could span a factor of ~ 50. 

Within this context it is interesting to compare the left-hand 
and right-hand panels of Fig. 5 where the limiting effect of restrict- 
ing the SED-fitting to a constant star-formation rate (CSF) model 
is explored. It can immediately be seen from the right-hand panel 
that if we adopt the same approach as Gonzalez et al. (201 1) and re- 
strict our SED-fitting to the CSF model, our stellar-mass estimates 
at z ~ 6.8 fall into excellent agreement with the AL, — Luv rela- 
tion they derived at z ~ 4. Moreover, it is also clear that restrict- 
ing the SED-fitting analysis to the CSF model significantly reduces 
(perhaps unrealistically) the scatter in the stellar-mass estimates at 
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Table 6. The best-fitting parameters returned by our SED-fitting analysis of the twenty-one objects in our final robust sample with detections at either 3.6/im or 
3.6 + 4.5^tm. The first column lists the objects IDs. The second column lists the best-fitting star-formation history (SFH) which is either: an instantaneous burst 
(Burst), constant star-formation rate (Const) or exponentially decaying star-formation rate (E). For the those objects where the best-fitting SFH is exponentially 
decaying (i.e. SFR <x e(~*/ T )), the proceeding number indicates the characteristic star-formation timescale in Gyrs (i.e. E0.2 => r = 0.2 Gyrs). Columns 3-8 
list the metallicity, age, reddening, stellar mass, star-formation rate and \ 2 of the best-fitting SED template. For those objects where the best-fitting SFH is a 
Burst, the SFR listed in column 7 is derived from the best-fitting model with a constant or exponentially decaying SFH. Note that the ages listed in columns 
4 & 10 refer to the total age of the system, rather than a luminosity-weighted age. Columns 9-12 list the x 2 , age, stellar mass and SFR of the best-fitting 
constant SFR model with a metallicity of O.2Z and A v = 0.0 (CSF model, see text for details). Based on this restricted set of SED parameters those objects 
highlighted with a f symbol in column 9 would have been rejected because their best-fitting SED templates have an unacceptably high x 2 - The final column 
lists an estimate of the SFR based on the UV luminosity of each object (M1500), with no correction for dust attenuation, which has been derived using the 
Madau, Pozzetti & Dickinson (1998) formula (corrected to a Chabrier IMF). 



ID SFH Z Age A v M* SFR xl xl Age M* SFR SFR UV 









(Myrs) 




(10 9 M ) 


(M yr x ) 






(Myrs) 


(10 9 Mq) 


(M yr x ) 


(M yr !) 


HUDF.1016 
HUDF.522 


Burst 
Burst 


0.02 
0.02 


10 
50 


1.1 

0.1 


Q+0.3 

1 O+0.9 


40.0±"'3 
15.01^ 


4.3 
4.5 


19.3t 
6.9 


575 
365 


q+0.5 

1 6+ 10 
1>o -0.4 


2 1+0.1 
za -o.i 
6 4+0-3 


3.8 
11.3 


HUDF.658 


Burst 


0.50 


65 


0.0 


1 Q+ 1 - 9 

ly -i.o 


c+3.3 
z -'-0.3 


1.4 


5.2 


725 


7 4+0.6 


5 1+0.1 
■'• 1 -0.4 


8.6 


HUDF.860 
HUDF.1173 


Burst 
Burst 


0.50 
0.02 


100 
55 


0.1 
0.0 




1 8+2-0 
1 1+ 2 ' 1 


1.8 
4.5 


17.8+ 
4.7 


645 
455 


i-5±C;S 
o.stgj 


, 4+0.3 

24+0.5 
Z - 4 -0.3 


5.4 
4.5 


HUDF09-2.1543 
HUDF09-2.2587 


Const 
Const 


0.02 
0.02 


725 
645 


0.1 
0.6 


2 1+3.1 
zl -1.3 


13.8± 4 7 4 6 


0.3 
3.3 


0.4 

9.5t 


575 
815 


1 5+ ' 9 
1 "'-0.3 

? 3+0.6 
z - J -0.9 


4 0+ ' 4 
4 1+ 01 

4a -0.3 


7.8 
7.2 


HUDF09-2.1660 


E1.0 


0.02 


725 


0.5 




10.0±|'f 


3.8 


8.7+ 


815 


2 ,+1.3 
za -0.5 


4 1+0-3 


7.2 


HUDF09-2.1584 


E0.2 


0.20 


725 


0.1 


16.5± 13 6 3 


3 n+0.5 


0.7 


7.6t 


645 


2 3+0.6 
z - J -0.9 


5 o+0.3 
J - J -0.3 


9.4 


ERS.7086 


Burst 


0.02 


25 


0.1 


0.4t° ;? 


4 9+4-1 


1.4 


1.7 


130 


4 +0A 


3 6+ 2 ' 3 


6.5 


ERS.6066 


E0.2 


1.00 


645 


0.1 


n 4+ 3 - 5 


4 4+0-4 


2.7 


25.4 1 " 


815 


24+0.6 
Z - 4 -0.9 


4 4+0-4 


7.2 


ERS.9100 


E0.2 


0.50 


405 


0.0 




2Q+1.3 


1.4 


3.2 


815 


1 4+0-4 


2 6 +0.3 
z -°-0.4 


4.5 


ERS.7225 
ERS.6438 


E0.05 
Burst 


0.20 
0.20 


325 
160 


0.0 
0.0 


4 °-2.3 


04+1-9 

3 6+ - 4 
J.O_ 4 


4.3 
6.7 


28.0t 
14.9+ 


815 
725 


1 9+0.5 
^-o.s 

2 ,+0.5 
za -0.8 


3 6 +0 ' 4 

4 5+0-3 


6.0 
7.2 


ERS.6263 


Burst 


0.02 


330 


0.0 


5 + 1 ' 1 
U.J_ 4 


<:+5.3 

z --'-o.8 


4.9 


5.3 


130 


4+°- 6 


4 3+4.4 


7.2 


ERS.7776 
ERS.5847 


E0.2 
E1.0 


0.02 
0.20 


325 
815 


0.0 
0.0 


1 1+ 11 
la -0.8 

3 o+4.9 


4.1^.8 


3.9 
2.7 


4.2 
2.9 


515 
815 


1 () +a6 

2 4+2.4 


9 9+O.5 
4 4+0-4 


5.4 
8.6 


ERS.3679 
ERS.7412 
ERS.6427 


Burst 
Burst 
E0.2 


0.02 
1.00 
0.02 


725 
10 
575 


0.2 
0.2 
0.0 


4 5+6-8 

4+0.8 

4 4+4-4 
^-1.6 


15 0+ ' 9 
'-'•"-ice 

11 4+71 

2 1+0.4 
Za -0.8 


5.5 
9.1 
1.2 


8.3 
12.2t 
6.6 


725 
255 
725 


4.3±i;i 
o.8±°;l 

2Q+0.5 


O9+0.3 
°- y -0.8 
39+I.I 

4 1+0-4 
4a -0.5 


16.4 
6.5 
7.2 


ERS.7376 


Const 


0.20 


455 


0.0 




3 

- , "-'-2.5 


0.6 


0.6 


455 


la -0.8 


3 5+0.9 
■ , - J -0.5 


6.5 



a given UV luminosity. Finally, as illustrated by the upper dotted 
line in the right-hand panel of Fig. 5, the brightest LBGs in our 
sample (i.e. M1500 < 19.0) are fully consistent with the expected 
M* — Luv relation for a galaxy which has been forming stars at 
a constant rate for ~ 600 Myrs. Importantly, at the mean redshift 
of the final robust sample (z = 6.8 ± 0.1), 600 Myrs represents 
~ 80% of the age of the Universe. Indeed, the primary cause of the 
clustering of objects around the M* — Ltjv relation corresponding 
to ~ 600 Myrs of constant star-formation is the requirement im- 
posed during the SED-fitting that objects must be younger than the 
age of the Universe. The underlying cause is simply that (with no 
dust reddening) the CSF models are bluer than the observed pho- 
tometry unless their age is close to the maximum allowable at this 
epoch. In summary, although the results shown in both panels of 
Fig. 5 are broadly compatible, it is clear that adopting a restricted 
set of SED templates may well provide a misleadingly low estimate 
of the true level of scatter in stellar mass at a given UV luminosity. 

Before moving on to consider the relationship between stellar 
mass and star-formation rate, it is worth remembering that one of 



the principal motivations for studying the M* — Ltjv relation at 
high redshift is to constrain the galaxy stellar-mass function (e.g. 
McLure et al. 2009). The results presented in Fig. 5 clearly illus- 
trate that in order to successfully constrain the stellar-mass func- 
tion at z > 6.5 is will be necessary to constrain the M* — Ltjv 
relation at UV luminosities substantially fainter than L*. Over the 
next three years, the new Cosmic Assembly Near-infrared Deep 
Extragalactic Legacy Survey (CANDELS; Co-P.I.s S. Faber & H. 
Ferguson; see Grogin et al. 2011 and Koekemoer et al. 2011) 
offers the prospect of significant progress. The deep portion of 
CANDELS will provide Y 105 , J125 &#i 6 u WFC3/IR imaging to 
uiab — 28(5ct) over an area of ~ 150 sq. arcmins in the GOODS 
N+S fields. The CANDELS programme should therefore provide a 
sample of > 200 robust z ~ 7 candidates in the magnitude range 
— 19 > M1500 > —20, all covered by the deep IRAC imaging 
available in GOODS N+S. A sample of this size should be suffi- 
cient to obtain robust constraints on the typical M*/Luv ratio at 
M1500 — — 19 by employing a stacking analysis to provide the nec- 
essary IRAC photometry. Obtaining constraints on the M* — Ltjv 
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relation at magnitudes as faint as ACoo — —18.5 (i.e. ~ 0.2L* 
at 2 ~ 7) will rely on stacking the final epoch2 WFC3/IR imaging 
of the HUDF into the forthcoming, ultra-deep, IRAC data being 
obtained as part of the Cycle 7 Spitzer warm mission GO-70145 
(P.I. Labbe). 

5.2 Star-formation rate versus stellar mass 

Over recent years it has become clear that studying the ratio of the 
current star-formation rate to the previously assembled stellar mass 
(SFR/M*), the so-called specific star-formation rate (sSFR), can 
provide a useful insight into the average star-formation history of a 
galaxy population. Studies of star-forming galaxies in the SDSS at 
z ~ 0.1 (Brinchmann et al. 2004), in the Extended Groth Strip at 
2 ~ 1.0 (Noeske et al. 2007) and in the GOODS fields at z ~ 1.0 
and z ~ 2.0 (Elbaz et al. 2007; Daddi et al. 2007) have consistently 
shown a correlation of the form SFR oc Af 09 (sSFR oc AC - 1 ), 
over a wide dynamic range in stellar mass. Moreover, these studies 
have shown that compared to a median value of sSFR ~ 2.5 Gyr -1 
at z = 2 (Daddi et al. 2007), the normalisation of the sSFR-M* 
relation has decreased by a factor of ~ 40 over the last 10 Gyrs. 

Indeed, there is some evidence that a sSFR of ~ 2.5 Gyr -1 
may represent the maximum sustainable star-formation rate for 
galaxies at 2 > 2. A recent study by Karim et al. (2011) found 
that high mass (M* > 10 M©) star-forming galaxies in the COS- 
MOS field consistently follow a steeper sSFR-M* relation than de- 
termined by previous studies (sSFR oc AC 0,4 ). However, by z ~ 3 
the results of Karim et al. (2011) suggest that the sSFR-M* rela- 
tion flattens at masses of M* < 10 Mq, becoming consistent with 
sSFR oc AC " 1 . Karim et al. (2011) interpret this flattening as a 
result of a natural limit to the sustainable sSFR of ~ 2.5 Gyr -1 , 
corresponding roughly to the inverse of the typical galaxy free-fall 
timescale. Within this context, it is interesting to note that Gonzalez 
et al. (2010) recently found that the median sSFR of a sample of 
twelve 2— drop candidates at 2 ~ 7 was also ~ 2.5 Gyr -1 . 

In Fig. 6 we plot SFR versus stellar mass for the twenty-one 
objects which have IRAC detections and therefore the most robust 
SFR and stellar-mass estimates. In the left-hand panel the stel- 
lar mass and dust corrected SFR estimates have been taken from 
the best-fitting SED templates drawn from the full range of star- 
formation histories, metallicities and dust reddening described in 
Section 3.3. In contrast, in the right-hand panel the stellar mass and 
SFR estimates have been taken from the best-fitting CSF model. 
Although the distribution of objects in the two panels is signifi- 
cantly different, both provide a consistent estimate for the typical 
sSFR. In the left-hand panel the median sSFR is 1.9 ± 0.8 Gyr -1 , 
while in the right-hand panel the median sSFR is 2.6 ± 0.4 Gyr - . 
Both estimates are clearly consistent with the typical sSFR value 
for 2 ~ 2 star-forming galaxies estimated by Daddi et al. (2007). 
To illustrate this point, in both panels of Fig. 6 the thick solid line 
is the best-fitting SFR-M* relation from Daddi et al. (2007) which 
corresponds to a sSFR of ~ 2.5 Gyr -1 . 

Consequently, taken at face value, our results provide addi- 
tional support to the conclusion that a direct proportionality be- 
tween SFR and stellar mass is still viable at 2 ~ 6.5, and that the 
corresponding sSFR of ~ 2.5 Gyr -1 may correspond to a physical 
limit on the maximum sustainable star-formation rate. However, it 
is clear from the left-hand panel that allowing a reasonable range 
of star-formation histories, metallicities and dust reddening leads 
to a large scatter in the SFR at a given stellar mass. Consequently, 
although the data shown in the left-hand panel are consistent with 
SFR and stellar mass being roughly proportional, they are also en- 



tirely consistent with star-formation and stellar mass being entirely 
unrelated. In contrast, the results shown in the right-hand panel sug- 
gest that SFR and stellar mass are well correlated, lying along a 
SFR— relation with a slope close to unity and a normalization 
consistent with a sSFR of ~ 2.5 Gyr -1 . 

However, it is worth noting that the apparently simple picture 
presented in the right-hand panel of Fig. 6 probably reflects limi- 
tations of relying on the CSF model, rather than offering genuine 
physical insight into high-redshift star-formation. The simple rea- 
son for this caution is that the agreement is largely inevitable when 
you only consider SEDs with constant star-formation and no red- 
dening. In this situation, each object is required to lie on a relation 
with slope of unity, with it's position on the SFR— M* plane simply 
determined by the best-fitting age. To illustrate this point we have 
plotted the expected SFR— M* relations for CSF models of various 
stellar population ages as the dotted lines in the right-hand panel 
of Fig. 6. This demonstrates that, provided the typical stellar popu- 
lation age lies in the range 200-600 Myrs, the resulting SFR— M* 
relation will automatically have a slope close to unity, and result in 
a typical sSFR consistent with ~ 2.5 Gyr -1 . 

In summary, although it is possible to constrain the typical 
sSFR of L* LBGs at 2 ~ 6.5, the limitations of the current sample 
do not allow meaningful constraints to be placed on the form of 
the SFR— M* relation. In order to resolve this issue it will be nec- 
essary to obtain much larger samples of 2 > 6 LBGs with stellar 
masses M* > 10 8 ' 5 M . Within this context, the new CANDELS 
WFC3/IR imaging data should prove decisive. The wide portion 
of CANDELS will proved J125 + ffieo imaging to a depth of 
m,AB — 27(5ct) over an area of ~ 0.2 sq. degrees, all of which 
is covered by deep IRAC imaging at 3.6 + 4.5/j.m {tuab — 26, 
5a) provided by the Spitzer Extended Deep Survey (SEDS; P.I. 
G. Fazio). The combination of CANDELS+SEDS should therefore 
allow the SFR— Af* relation at 2 ~ 6.5 to be investigated using 
a sample of >250 LBGs with reliable stellar-mass estimates of 
Af* > 10 8 ' 5 M Q . 

5.3 The effect of nebular emission 

Recent work has suggested that nebular continuum and line emis- 
sion might contribute to the observed SEDs of 2 ~ 6 — 7 galaxies 
(e.g. Ono et al. 2010; Schaerer et al. 2010). As discussed by Robert- 
son et al. (2010), galaxies with strong UV continua can typically be 
fit using pure stellar populations with ages of a few hundred Myr, or 
by much younger populations (< few Myr) with significant nebular 
contributions and an implied low escape fraction (/ es c) of Lyman 
continuum photons. Such nebular solutions can yield much lower 
stellar masses than those in the purely stellar case (Ono et al. 2010). 

In order to quantify this degeneracy and its possible effect on 
our derived physical properties, we have examined in more detail 
the sub-sample of twenty one galaxies detected in the 3.6 fim IRAC 
band (nine of which are also detected at 4.5/im). For these objects it 
is possible to investigate whether the IRAC detections can provide 
a valuable discriminant between the nebular and stellar solutions, 
given the location of prominent nebular lines, such as H/3 and [OIII] 
5007 A, at the redshifts of interest. 

5.3.1 Nebular emission methodology 

As before we use the Bruzual & Chariot (2003) models to gener- 
ate a set of spectral templates based on a Chabrier IMF. For the 
models presented here, we use a representative exponentially de- 
caying star-formation history (r ~ 0.4 Gyr) and, motivated by the 
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Figure 6. Star-formation rate versus stellar mass for the twenty-one objects in our final robust sample with IRAC detections at either 3.6/jm (grey) or 3.6+4.5/jm 
(black). In the left-hand panel the star-formation rates and stellar masses have been measured from the best-fitting SED template drawn from the full range 
of star-formation histories, metallicities and reddening described in Section 3.3. In the right-hand panel the star-formation rates and stellar masses have been 
estimated from the best-fitting CSF model (see text for details). The la errors on both parameters have been calculated by determining the A\ 2 = 1 interval, 
after marginalization over all other free parameters. The solid line in both panels is the SFR— M* relation derived by Daddi et al. (2007) for star-forming 
galaxies at z ~ 2 and corresponds to a sSFR of ~ 2.5 Gyr - 1 . In the right-hand panel the dotted lines illustrate how the SFR— M* relation for a galaxy with 
a constant star-formation rate and zero reddening varies as a function of stellar population age. 
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Figure 7. An example of our SED-fitting procedure using models which 
incorporate nebular continuum and line emission. In this illustration the 
pure stellar continuum model (/ CS c = 1) is shown as the red (thick) line 
and the template featuring nebular emission (/ eS c = 0.2) is shown as the 
blue (thin) line. For clarity the galaxy template featuring nebular emission 
has been displaced by 0. 1 magnitudes (both models have a metallicity of 
Z = O.2Z ). 



L(Hf3) = 4.78 x HT 13 (1 - f csc )N lyc (6) 

Other HI line intensities follow from ratios predicted by standard 
recombination theory (Osterbrock & Ferland 2006). Lines from 
common metallic species are included using relative intensities 
given by Anders et al. (2003) assuming the gas phase metallicity is 
either Zq or 0.2 Zq. We use the method of Brown et al. (1970) to 
calculate the strength of bound-free and free-free continuum emis- 
sion, and use results from Osterbrock & Ferland (2006) for the two 
photon emission from H. 

When fitting the SED models we consider two fixed values of 
the escape fraction, / csc = 0.2 (stellar and nebular emission) and 
/esc = 1 (purely stellar emission). The value of / csc = 0.2 is moti- 
vated by direct observations of the Lyman continuum in galaxies at 
z ~ 3 (Shapley et al. 2006) and typical values of / csc required for 
star-forming galaxies to maintain reionization at z ~ 7 (Robertson 
et al. 2010). As we are primarily interested in how nebular emission 
might alter the inferred stellar mass and age, we do not include the 
possible effects of Lya emission or reddening. An example SED fit 
featuring nebular emission is shown in Fig. 7 and the 0.2 Zq nebu- 
lar SED fits for all twenty one objects can be found in Appendix B. 



fact that high-redshift galaxies often exhibit low metallicities (e.g. 
Finkelstein et al. 2011), we consider both solar (Z = Zq) and fifth 
solar models (Z = 0.2 Zq). In addition to exponentially decaying 
star-formation histories, we also investigated models with constant 
star formation, but found that these did not significantly alter our re- 
sults. The contribution of the nebular continuum and line emission 
is computed in the manner of Robertson et al. (2010), providing 
nebular emission models similar to those calculated by Ono et al. 
(2010). The strength of the nebular emission is tied to the number 
of ionizing photons per second (A?Lyc), calculated from the stellar 
population model, via the H/3 luminosity (in erg s _1 ): 



5.3.2 Nebular emission results 

For the models with sub-solar metallicity it is found that, compared 
to the purely stellar templates (/ csc = 1), the templates which in- 
clude nebular emission (/ csc = 0.2) provide a better fit to 14 of 
the 21 objects, and all nine of the objects with detections in both 
IRAC bands. The data favour the nebular emission models because 
of the predominantly blue [3.6] - [4.5] fim observed colours that are 
easily reproduced by including rest-frame optical line emission. In 
contrast, the shape of the purely stellar model SEDs redward of the 
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Figure 8. A histogram of the ratio of the best-fitting stellar masses returned 
by fitting Z = 0.2 Z Q SED models with (M n ) and without (M*) nebular 
continuum and line emission for the twenty one objects with 3.6/^m IRAC 
detections. For 16 out of 21 objects the stellar-mass estimates differ by less 
than a factor of 2.5. However, it can be seen that there are two objects for 
which the nebular fits return stellar masses which are an order of magnitude 
smaller than the stellar masses returned by fits without nebular emission. 

Balmer/4000A break cannot easily accommodate these blue rest- 
frame optical colours. 

For the majority of this sub-sample (16 out of 21) the best- 
fitting ages are still > 100 Myr with stellar masses reduced by less 
than a factor of 2.5 (see Fig. 8). In contrast, the best-fitting mod- 
els to the remaining five galaxies have significantly lower stellar 
masses (by more than a factor of ten in two cases). Interestingly, 
all five of these objects are only detected in the 3.6 /im band and 
have best-fitting ages which are perhaps unphysical (median age 9 
Myrs). Using solar metallicities, we find eight galaxies with signif- 
icantly altered parameters, including all five found for models with 
Z = O.2Z . 

In summary, with the current precision of the HST and IRAC 
photometry, we are unable to draw firm conclusions about the pos- 
sible presence of nebular emission in these sources. However, given 
that the inferred stellar masses from the models including nebular 
emission are generally similar to those inferred from purely stel- 
lar models, we conclude that the results presented in Sections 5.1 
and 5.2 appear to be robust to the inclusion of / csc = 0.2 nebular 
emission, especially for those objects with detection in both IRAC 
bands. 



6 COMPARISON WITH PREVIOUS STUDIES 

As discussed in the introduction, the availability of the various new 
WFC3/IR datasets has led to a proliferation of papers focused on 
z > 6 LBGs. With authors each applying their individual candidate 
selection procedures, and in many cases using their own indepen- 
dent reductions of the publicly available data, it is difficult to obtain 
a clear overview of the subject, and to identify whether different 
studies are in good agreement, or not. Consequently, in this section 
we compare our final robust sample with those derived elsewhere 
in the literature (on a field-by-field basis), highlighting the objects 
we have in common, and investigating the properties of previously 



published high-redshift candidates which are not included in our 
final robust sample. 

Given that one of the primary motivations for this study was 
to derive a sample of high-redshift candidates which is as robust 
as possible, and that previous samples of WFC3/IR high-redshift 
candidates were selected for a variety of purposes, it is not the 
case that we regard any object not included in our final sample 
as a low-redshift interloper. In fact, as the proceeding discussion 
will demonstrate, each of the candidates from the literature sam- 
ples falls into one of four categories. The first category consists 
of objects which are in common with our final sample of seventy 
z > 6 LBGs and we therefore regard as being robust. The second 
category consists of objects which were not included in our final 
sample (because they failed to meet one or more of our adopted 
criteria), but which nevertheless our analysis suggests are likely to 
be at high redshift. The third category consists of objects which 
our analysis suggests are likely to be a low redshift, but do have 
an acceptable (albeit lower probability) solution at high redshift. 
The fourth category consists of those objects which our analysis 
suggests are very unlikely to be at high redshift. Throughout the 
discussion in this Section we have attempted to make it as clear 
as possible which category each of the candidates falls into. Fi- 
nally, it should be noted that where a research group has published 
a number of studies of a particular survey field, we only discuss the 
results from the most recent study, under the assumption that they 
supersede any previous work. 



6.1 HUDF 

6.1.1 McLure et al. 2010 

In McLure et al. (2010) we published our initial analysis of the 
HUDF WFC3/IR dataset, providing a list of N=49 high-redshift 
candidates with z p hot > 5.9. The candidate selection procedure 
employed in McLure et al. (2010) was broadly similar to that 
adopted here, with the most noteworthy difference between the two 
analyses being that in this work we have directly employed decon- 
fused IRAC photometry in the candidate selection procedure. 

Of the N=31 objects identified in the final HUDF sample listed 
in Table 2, N=28 are in common with the sample derived in McLure 
et al. (2010), demonstrating an excellent level of agreement be- 
tween the two studies. However, there are N=21 candidates pub- 
lished in McLure et al. (2010) which do not feature in the final 
robust sample derived here. The reason behind this is that primary 
aim of McLure et al. (2010) was to provide an estimate of the z = 7 
and z = 8 galaxy luminosity functions. Consequently, the McLure 
et al. (2010) sample was designed to be as complete as possible, 
and therefore contained all potential z > 5.9 LBGs revealed by the 
SED-fitting analysis, irrespective of whether or not they also dis- 
played an acceptable low-redshift solution^. In contrast, the princi- 
pal aim of this study is to derive a sample of z > 6 LBG candidates 
which is as robust as possible, which means, in effect, requiring that 
any alternative low-redshift solutions can be statistically excluded. 
Although all N=21 of the additional candidates listed in McLure et 
al. (2010) also feature in the initial catalogues derived here, all of 
them were excluded from our final robust HUDF sample because 
the best-fitting alternative solutions at low-redshift could not be ex- 
cluded at the Ax 2 > 4 (95%) level. 



5 Note that the alternative low-redshift solutions were also listed by 
McLure et al. (2010). 
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6.1.2 Bouwens et al. 2011 

Based on their analysis of the HUDF dataset, Bouwens et al. (201 1) 
list a total of N=31 robust high-redshift candidates, which are a 
mixture of z — drops and Y — drops. Of these N=31 candidates, eigh- 
teen also featured in our final robust sample of HUDF candidates 
listed in Table 2. However, it is clearly of interest to investigate why 
the remaining thirteen objects identified by Bouwens et al. (2011) 
do not feature in our final robust sample. 

Firstly, we should note that six of the thirteen additional 
objects (UDFz-38537518, UDFy-37588003, UDFy-33446598, UDFy- 
39347255, UDFy-40338026 & UDFy-42406550) are simply too faint to 
make it into our final robust sample. None of these six objects is 
bright enough (in a 0.6"diameter aperture) to provide a > 5cr de- 
tection in any of the WFC3/IR bands, and all six are fainter than 
any of our robust HUDF candidates. Consequently, this leaves a to- 
tal of seven high-redshift candidates listed by Bouwens et al. (201 1) 
which could, in principle, also feature in our final robust sample. 

Our SED-fitting analysis suggests that two of the additional 
objects (UDFz-44746449 & UDFy-43086276) are likely to be at high- 
redshift (z p hot = 8.1 and z v hot = 8.3 respectively), but just failed 
to make our final robust sample because the competing low-redshift 
solutions could not be ruled-out at > 95% confidence. A further 
two additional objects (UDFz-42567314 & UDFz-42247087) also have 
primary photometric redshift solutions at z p hot > 6.0, but were 
subsequently rejected because they were either too close to the 
WFC3/IR array edge (UDFz-425673l43), or were deemed to have 
unreliable photometry due to contamination from a nearby, bright, 
low-redshift galaxy (UDFz-42247087). Of the remaining three ob- 
jects, one (UDFy-37796001) does have an acceptable solution at 
Zphot > 8, but was rejected because our analysis suggests that the 
alternative solution at z v hot — 2 is marginally preferred. The other 
two (UDFz-37296175 & UDFy-37636015) were rejected because our 
SED fitting analysis returned a primary photometric redshift solu- 
tions at z phot — 5. 

Finally, it can be seen from Table 2 that our final robust sam- 
ple contains thirteen objects which are not featured in the Bouwens 
et al. (201 1) robust candidate list. However, the noteworthy feature 
of these objects is that the vast majority (11/13) are at z p hot < 
6.5, whereas the Bouwens et al., colour-colour, selection criteria 
are tuned to select objects at z p hot > 6.5. The two exceptions 
(HUDF.1730 & HUDF.2701) have been identified by several dif- 
ferent studies (see Table 2 for details) and one (HUDF.2701) does 
feature in the Bouwens et al (2011) list of potential high-redshift 
candidates. 



6. 1.3 Finkelstein et al. 2010 

In their analysis of the WFC3/IR HUDF dataset, Finkelstein et al. 
(2010) used a similar template-fitting technique to that employed in 
both McLure et al. (2010) and this work, and used each candidate's 
photometric redshift probability density function in the construc- 
tion of their final list of N=31 candidates at 6.3 < Zphot < 8.6. 
As part of their analysis, Finkelstein et al. (2010) conducted a de- 
tailed comparison between their final list of high-redshift candi- 
dates and the McLure et al. (2010) sample, finding a good level of 
agreement between the two studies. 

As might be expected, the overall agreement between the anal- 
ysis of Finkelstein et al. (2010) and the final robust HUDF sample 



reported as ID=1 144 in McLure et al. (2010) 



derived here is still good. In the redshift range covered by both 
studies, our final robust HUDF sample consists of N=22 candidates 
at Zphot > 6.3, eighteen of which are in common with Finkel- 
stein et al. (2010). The four additional candidates which feature in 
our final robust sample are: HUDF.2281, HUDF.2324, HUDF.2672 & 
HUDF.2664 (see Appendix B for plots of the SED fits). 

Of the N=31 candidates in the Finkelstein et al. (2010) sample, 
N=19 also feature in the final HUDF sample derived here. However, 
this still leaves a total of twelve candidates from Finkelstein et al. 
(2010) which do not feature in our final sample. All twelve of these 
additional candidates do feature in our original HUDF catalogues, 
but were excluded from the final robust sample for a number of dif- 
ferent reasons. One object (FID 3022) was excluded from our sam- 
ple because it is too faint (J125 > 29) to provide a robust high- 
redshift solution, and a further four objects (FIDs 640, 1818, 2013 & 
2432) were excluded because they were judged to have photometry 
which was potentially contaminated by bright, nearby, low-redshift 
galaxies. For the remaining seven objects (FIDs 200, 213, 567, 653, 
1110, 1566, & 2055) our SED fitting analysis does indicate that the 
primary photometric redshift solution is at z p hot > 6.3. However, 
all seven objects were excluded from the final robust sample be- 
cause our analysis suggested that the alternative low-redshift solu- 
tion could not be ruled-out with > 95% confidence. 

6.1.4 Yan et al. 2010 

In their analysis of the HUDF, Yan et al. (2010) used z— drop and 
Y— drop criteria to identify a sample of N=35 high-redshift can- 
didates at z ~ 7 and 2 ~ 8. Excluding a likely transient, Yan 
et al. (2010) list a total of twenty drop candidates, fourteen of 
which are in common with our final robust HUDF sample. Of the 
six z— drops listed by Yan et al. (2010) which don't make it into 
our final robust HUDF sample, two (A046 & A056) were excluded 
because their alternative low-redshift solutions could not be ruled- 
out with > 95% confidence, one (A017) was excluded because its 
photometry was contaminated by a bright, low-redshift, galaxy and 
one (A008) was rejected because it lies too close to the array edge. 
The final two z-drops (A055 & A062) listed by Yan et al. (2010) do 
not feature in any of our catalogues and do not appear to be robust 
objects based on our reduction of the epoch 1 HUDF dataset. 

Yan et al. (2010) list a total of fifteen Y — drop candidates in 
the HUDF. Of these fifteen candidates, only two (B092 & B115) 
make it through to our final robust sample. Of the thirteen Y — drops 
listed by Yan et al. (2010) which do not feature in our final sam- 
ple, our analysis suggests that five (B041, B088, B114, B117 & SB27) 
do have acceptable high-redshift photometric redshift solutions, but 
were excluded because they all have alternative low-redshift solu- 
tions which cannot be ruled-out at the >95% confidence level. Two 
further objects (B087 & B094) also feature in our original catalogues 
but, based on our 0.6"diameter photometry, are not Y — drops and 
have primary photometric redshift solutions at z p hot < 7- The re- 
maining six candidates (SB30, SD02, SD05, SD15, SD24 & SD52) do 
not appear as robust objects in our reduction of the epoch 1 HUDF 
dataset. Finally, we note that Yan et al. (2010) also identify a sam- 
ple of twenty three J— drops in the HUDF, none of which feature 
in our final robust HUDF sample. 

6.1.5 Wilkins et al. (2010) 

Wilkins et al. (2010) identify a total of eleven z— drop candidates 
in the HUDF, nine of which also feature in our final robust sam- 
ple. Of the two additional candidates listed by Wilkins et al., our 



© 0000 RAS, MNRAS 000, 000-000 



16 R.J. McLure et al. 




Figure 9. Stacked postage-stamp images of three z ~ 7 candidates (ERS.z.80252, ERS.z.47667 & ERS.z.20851) published by Wilkins et al. (2010). All 
three of these objects were rejected as high-redshift candidates by our SED fitting analysis, principally due to low signal-to-noise detections in the blue optical 
bands. This figure demonstrates that although there is a drop in flux between the 2850 and Yrj98 filters, these objects are also clearly detected in a stack of the 
-B435 + Veoe + *775 GOODS imaging, suggesting that they are unlikely to be at z ~ 7. 



analysis suggests that one (HUDF.z.6497) does have an acceptable 
solution at z p hot > 6, but was was excluded because the primary 
photometric redshift solution lies at z v hot = 3.5. The other object 
(HUDF.z.6433 6 ) was rejected because it lies close to the array edge 
and was therefore deemed to have unreliable photometry. 

6.7.6 Lorenzoni et al. (2011) 

Based on their analysis of the HUDF dataset, Lorenzoni et al. 
(2011) identify a sample of six Y— drop candidates. Of these six 
candidates, three (HUDF.YD1, HUDF.YD3 & HUDF.YD4) make it 
into our final robust HUDF sample. Of the remaining three can- 
didates, our SED fitting analysis suggests that two (HUDF.YD2 & 
HUDF.YD8) have an acceptable z ~ 8 photometric redshift so- 
lution, but were excluded from our final robust sample because 
they both have an alternative low-redshift solution which cannot 
be securely ruled out (i.e. A\ 2 < 4). The remaining candidate 
(HUDF.YD9) does not appear as a robust object in any of our cata- 
logues. 

6.2 ERS 

6.2.1 Bouwens et al. (201 1 ) 

The robust ERS sample derived by Bouwens et al. (2011) con- 
sists of N=19 objects in total, thirteen z~ drops at 2 ~ 7 and six 
Y— drops at z ~ 8. Of the thirteen z— drops listed by Bouwens et 
al. (201 1), only five appear in our final robust sample (see Table 3). 
Of the remaining eight additional z— drops listed by Bouwens et al. 
(2011), one object (ERSz-2352941047) is too faint (J125 > 27.5) 
to produce a robust high-redshift solution based on our criteria, 
leaving seven additional z— drops to account for. Of these, three 
(ERSz-2 150242362, ERSz-2225141173 & ERSz-2354442550 ) have sta- 
tistically acceptable photometric redshift solutions at z p hot > 6.5, 
but were excluded from the final robust sample because it was not 
possible to rule-out the alternative low-redshift solutions at > 95% 
confidence. One further object (ERSz-2150943417) was rejected be- 
cause based on our photometry it wasn't possible to obtain a sta- 
tistically acceptable solution at high-redshift. Of the final three ob- 
jects, two (ERSz-21 11644168 & ERSz-2432842478) have acceptable 
solutions at z p hot > 6 but were excluded because our primary 
photometric redshift solution lies at z p hot < 2. The final object 
(ERSz-2056344288) does not have an acceptable high-redshift solu- 
tion based on our 0.6"diameter aperture photometry. 

Of the six Y— drops listed by Bouwens et al. (2011), two 
(ERSY-2354441327 & ERSY-2029843519) make it into are final robust 
sample. Of the four additional Y — drops listed by Bouwens et al. 
(201 1), one object (ERSY-2377942344) is too faint in a 0.6" diameter 



aperture ( J125 > 27.5) to produce a robust high-redshift solution 
based on our criteria, leaving three additional Y — drops to be ac- 
counted for. Of these three objects, two (ERSY-2399642019 & ERSY- 
2251641574) have acceptable primary photometric redshift solutions 
at z > 7.5, but were excluded because the alternative low-redshift 
solutions could not be rule-out. The final object, ERSY-23Q6143O40 
was rejected because our primary photometric redshift solution is 
sit z p hot ^ 2. 

6.2.2 Lorenzoni et al. (2011) 

Lorenzoni et al. (201 1) identify a total of nine Y — drop candidates 
in the ERS field (five of which, marked with *, are described as 
"more marginal candidates") . Of these nine candidates, only two 
(ERS.YD1 & ERS.YD2*) make it into our final robust sample. Of 
the remaining seven candidates, our analysis suggests that three 
(ERS.YD5*, ERS.YD6 & ERS.YD9*) have an acceptable solution at 
z p hot > 7, but were rejected because the alternative low-redshift 
solution could not be ruled out at > 95% confidence. A further 
two candidates (ERS.YD7* & ERS.YD8*) were excluded because 
our primary photometric redshift solution lies at z p hot — 2. The 
remaining two objects (ERS.YD3 & ERS.YD4) do not appear as ro- 
bust objects in any of our catalogues. 

6.2.3 Wilkins et al. (2010) 

Based on their analysis of the ERS field, Wilkins et al. (2010) iden- 
tify a sample of eleven z— drop candidates, six of which also feature 
in our final robust sample. Of the five additional candidates listed 
by Wilkins et al., one object (ERS.z.26813) does have an accept- 
able primary photometric redshift solution at z p hot = 6.6, but was 
excluded from our final sample because it has an equally accept- 
able solution at z v y tot — 1.5. A further object (ERS.z.70546) was 
rejected because it was not possible to obtain an acceptable high- 
redshift SED fit. The three remaining candidates listed by Wilkins 
et al. (ERS.z.80252, ERS.z.47667 & ERS.z.20851) were rejected as 
low-redshift interlopers by our SED-fitting analysis due to the pres- 
ence of consistent, low-level, detections in the bluer optical bands. 
To illustrate this point we have stacked the ACS+WFC3/IR data for 
these three objects and show the resulting postage-stamp images in 
Fig. 9. It can clearly be seen that although there is a drop in flux 
between the zsso and Y"o98 filters, the significant detection of flux 
in the stack of the B435 + V606 + *775 images suggests these objects 
are unlikely to be at z ~ 7. 

7 This object was highlighted by Bouwens et al. (201 1) as being potentially 
at low redshift. 
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6.3 HUDF09-2 

6. 3. 1 Bouwens etal. (201 1 ) 

Bouwens et al. (2011) lists a total of N=35 robust high-redshift 
candidates in the HUDF09-2 field, consisting of eighteen z-drops 
and seventeen Y-drops. Only seven of these thirty-five candidates 
appear in our final robust sample (including HUDF09-2.799 which 
requires a contribution from Lya line emission), which clearly re- 
quires some explanation. The principal reason for this apparent dis- 
crepancy is that the Bouwens et al. (2011) sample contains many 
fainter objects than our final robust sample. Indeed, of the thirty- 
five candidates listed by Bouwens et al., seventeen are fainter (in 
our 0.6"diameter aperture photometry) than the faintest member of 
our final robust sample. Therefore, based on the data utilised in this 
study, and our criteria for isolating robust candidates, it is likely 
that these seventeen objects are simply too faint to make it into our 
final robust sample. 

However, even accounting for the difference in selection 
depth, there are still eleven robust candidates identified by 
Bouwens et al. (2011) which should, in principle, also appear in 
our final robust sample. All eleven of these candidates do feature 
in our HUDF09-2 catalogues, but were excluded from the final ro- 
bust sample for a number of reasons. Three of the additional candi- 
dates (UDF092z-008 11320, UDF092z-07091160 & UDF092y-07090218) 
have acceptable high-redshift photometric redshift solutions, and 
were close to making our final robust candidate list. However, 
for these candidates, the difference in x 2 between the primary 
photometric redshift solution and the alternative low-redshift so- 
lution (Ax 2 — 3) did not quite match our adopted criterion of 
Ax 2 > 4. Of the remaining eight additional candidates listed by 
Bouwens et al. (2011), five (UDF092y-02731564, UDF092z-09770485, 
UDF092z-09151531, UDF092y-06321217 & UDF092y-06391247) were 
rejected because our primary photometric redshift solutions lie 
in the redshift interval 4.9 < z p h ot < 5.9. The remaining 
three additional candidates (UDF092y-04242094, UDF092y-09611126 
& UDF092y-09661163) were rejected because our analysis suggests 
that their primary photometric redshift solutions are at z p hot — 2.1. 

Finally, we should note that two of the candidates which ap- 
pear in our final robust sample (HUDF09-2.2455 & HUDF09-2.2814) 
feature in the Bouwens et al. (201 1) list of potential, but non-robust, 
high-redshift candidates. Moreover, our final robust sample fea- 
tures eight candidates which do not appear in any of the Bouwens 
et al. (2011) lists, although 7/8 of these additional candidates have 
Zphot < 6.3, where the z— drop criteria applied by Bouwens et al. 
is less sensitive. 



6.3.2 Wilkins et al. (2010) 

Wilkins et al. (2010) list a total of fifteen z— drop candidates in the 
HUDF09-2 field, of which six also appear in our final robust sam- 
ple. All nine of the additional candidates listed by Wilkins et al. 
(2010) feature in our initial HUDF09-2 catalogues, but were ex- 
cluded from the final sample for a number of different reasons. 
One candidate (P34.Z.3996) does have a valid high-redshift solu- 
tion at Zph t = 6.55, and was close to making it into the fi- 
nal robust sample, but was rejected because the A^ 2 between the 
high-redshift and low-redshift alternative solution was too small 



8 Bouwens et al. (201 1) exploit deep F814W imaging which partially cov- 
ers the HUDF09-2 field and, in some cases, will allow the selection of 
fainter high-redshift candidates. 




log(A/7tm) 




log(X/^m) 




log(X///m) 

Figure 10. The results of our SED-fitting analysis for three z ~ 7 candi- 
dates in HUDF09-2 (P34) from Wilkins et al. (2010). Each of these objects 
was rejected from our final robust sample because our SED-fitting analysis 
does not return a plausible high-redshift solution. In each panel the solid 
(red) line shows the best-fitting galaxy template at low redshift while the 
dotted (blue) line shows the alternative high-redshift solution (note: there 
is no high-redshift solution for P34.Z.703). The best-fitting photometric 



redshifts for these candidates are 



1.4, 1.6&1.6 for P34.Z.703, 



z phot 

P34.Z.4288 and P34.Z.5016 respectively. In each panel the inset showing 
X 2 as a function of redshift demonstrates that any solution at z ~ 7 has a 
very low probability. 



(Ax 2 = 3.4). Of the remaining eight candidates, two candidates 
(P34.Z.703 & P34.Z.2428) were excluded because their primary pho- 
tometric redshift solutions are in the interval 4.5 < z v hot < 5.5. 
The remaining six candidates (P34.Z.2397, P34.Z.3053, P34.Z.4288, 
P34.Z.4501, P34.Z.5016 & P34.Z.3990) were rejected because their pri- 
mary photometric redshift solutions were all at z p hot < 4.5. As an 
illustration, in Fig. 10 we show our SED fits for three of the 2— drop 
candidates from Wilkins et al. (2010). In each case our analysis sug- 
gests that there is a very low probability of the candidate being at 
z ~ 7. 
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6. 3. 3 Lorenzoni et al. (2011) 

Lorenzoni et al. (2011) list a total of seven Y— drop candidates in 
the HUDF09-2 field (two highlighted with a * as more marginal 
candidates). None of these seven candidates make it into our fi- 
nal robust sample. Four of the seven Lorenzoni et al. Y — drop 
candidates (P34.YD1, P34.YD2, P34.YD3 & P34.YD6*) feature in our 
original HUDF09-2 catalogues, and do have acceptable photomet- 
ric redshift solutions at 7.0 < z p hot < 8.5. However, all four 
objects have comparable/preferred solutions in the redshift range 
1.3 < Zphot < 2.1 and were therefore rejected. Of the remaining 
objects, one (P34.YD7*) is not a robust object in our reduction of 
the HUDF09-2 dataset. The final two candidates were excluded be- 
cause they were deemed to have unreliable photometry, either due 
to being within the wings of a bright star (P34.YD4) or for being too 
close to the noisy array edge (P34.YD5). 



7 CONCLUSIONS 

We have presented the results of a study designed to identify robust 
high-redshift (z > 6) galaxies using the available multi-wavelength 
(ACS+WFC3/IR+IRAC) imaging covering the HUDF, HUDF09- 
2 and ERS fields. By exploiting the advantages provided by our 
SED-fitting analysis, and by incorporating deconfused IRAC pho- 
tometry directly into our candidate selection procedure, we have 
produced a robust sample of seventy high-redshift candidates at 
redshifts 6.0 < z p hot < 8.7. Based on this final robust sample 
we have investigated the distribution of UV spectral slopes, stellar 
masses and star-formation rates. Finally, we have presented the re- 
sults of a detailed comparison between our final robust sample and 
previous samples of high-redshift candidates derived from the same 
datasets using different selection criteria. Our main conclusions can 
be summarised as follows: 

(i) Based on our final sample of seventy robust high-redshift 
candidates, and employing a variance weighted mean to account 
for the wide range in individual uncertainties, we find that the typ- 
ical value of the UV spectral slope is (/3) = —2.05 ± 0.09. Conse- 
quently, in contrast to some previous studies, we find no evidence 
that LBGs in the redshift interval 6.0 < z p ) lo t < 8.7 display 
UV spectral slopes which are significantly bluer than seen in star- 
forming galaxies at lower redshifts. 

(ii) Using the stellar-mass estimates provided by our SED-fitting 
we have investigated the relationship between stellar mass and UV 
luminosity at z ~ 7, finding a best-fitting relationship of the form: 
M* oc Lijoo ' 6 . We note that our best-fitting M* - L vv relation 
is consistent with, although somewhat steeper than, the M* — Luv 
relations derived by previous authors based on large samples of 
LBGs at z ~ 4. 

(iii) Focusing on a sub-sample of high-redshift candidates 
with reliable IRAC photometry, we find that L* LBGs at 
z ~ 6.5 (i.e. Misoo — —20.2) have a median stellar mass of 
Mi, = (2.1 ± 1.1) x 1O 9 M . However, by employing SED 
templates featuring a range of metallicities, star-formation histo- 
ries and reddening we find that the scatter in stellar mass at a given 
UV luminosity could span a factor of ~ 50. We caution that de- 
riving stellar-mass estimates solely from models with a constant 
star-formation rate may provide a misleadingly low estimate of the 
real range in stellar mass at a given UV luminosity. 

(iv) Based on the same sub-sample of twenty one objects with 
IRAC detections, we find that the median specific star-formation 
rate (sSFR) of L* LBGs at z ~ 6.5 is 1.9 ± 0.8 Gyr -1 , a value 



which is consistent with the sSFR ~ 2.5 Gyr -1 observed in star- 
forming galaxies at z ~ 2. However, our SED fitting indicates that 
the range in the star-formation rate at a given stellar mass is po- 
tentially large and that adopting constant star-formation rate mod- 
els with zero reddening may provide an underestimate of the real 
level of scatter. Moreover, we find that SED fitting based solely on 
constant star-formation rate models (with zero reddening) will in- 
evitably tend to produce a typical sSFR close to ~ 2.5 Gyr -1 for 
high-redshift galaxies with ages in the range 200 — 600 Myrs. 

(v) Using the sub-sample of twenty one objects with IRAC de- 
tections we have also investigated the impact on our stellar-mass es- 
timates of including the effects of nebular continuum and line emis- 
sion in our galaxy SED templates. Based on SED templates with 
sub-solar metallicity (Z = 0.2 Z Q ) we find that nebular (/ csc = 0.2) 
models are typically capable of providing a better fit to the predom- 
inantly blue [3.6] - [4.5] pirn colours. However, in the majority of 
cases (16 out of 21) the best-fitting stellar masses returned by the 
nebular fits are less than a factor of 2.5 lower than those returned 
by the stellar continuum fits. Those objects for which the nebular 
fits return significantly lower stellar masses (5 out of 21) are found 
to have a somewhat unphysical median age of 9 Myrs. 

(vi) A detailed comparison between our final robust sample 
and previous high-redshift samples derived using different selec- 
tion/analysis techniques produces mixed results. When confined to 
the brightest objects, and the best data, the overlap between our 
final robust sample and samples derived using traditional "drop 
out" criteria is reasonably good. However, at fainter magnitudes 
our analysis suggests that some literature samples are heavily con- 
taminated (in some cases at the > 50% level) by objects which are 
very likely at lower redshifts. 
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APPENDIX A: IRAC DECONFUSION ALGORITHM 

As previously discussed in Section 2.3, deep Spitzer IRAC imag- 
ing data is available for all three of the fields analysed in this 
paper (3.6/im+4.5/im for the HUDF+ERS fields, and 3.6/im for 
HUDF09-2) with total integration times ranging from 23 to 46 
hours. Fully exploiting the information provided by the IRAC 
imaging is vital for analysing the high-redshift galaxy population 
for two fundamental reasons. Firstly, given that even the Hiqo 
imaging data is sampling rest-frame wavelengths of A < 2200A at 
z > 6, the information long- ward of the 4000A break provided by 
the IRAC imaging is crucial for providing constraints on the stellar 
mass and the contribution of any older stellar population (see Sec- 
tion 5). Secondly, inclusion of the IRAC photometry in the SED 
fitting process is very effective at breaking photometric redshift de- 
generacies between genuine z > 6 candidates, and interlopers at 
both 2 ~ 2 and z ~ 5 (see SED plots in Appendix B). 

Although the availability of ultra-deep IRAC imaging is po- 
tentially hugely beneficial, because the IRAC data is heavily con- 
fused, obtaining accurate flux measurements for faint, high-redshift 
galaxies is technically challenging. Several approaches to solving 
this problem have been presented in the literature, all of which rely 
on using a higher resolution image (ideally as close as possible in 
wavelength) as prior information to deconfuse the lower-resolution 
IRAC data. One approach, which has been recently applied to the 
ultra-deep IRAC data in the GOODS-N field, uses a modified ver- 
sion of the CLEAN algorithm (Hogbom 1974), traditionally used 
in radio astronomy, together with model templates extracted from 
high-resolution, ground-based, K— band imaging to estimate IRAC 
fluxes via an iterative subtraction scheme (Wang et al. 2010). An al- 
ternative approach relies on building two-dimensional, axisymmet- 
ric models of each the galaxy based on the available high-resolution 
imaging (e.g. Labbe et al. 2006). The two-dimensional models are 
then convolved to the spatial resolution of the IRAC imaging and 
a x 2 minimization procedure is employed to determine the indi- 
vidual galaxy fluxes which best reproduce the observed IRAC im- 
age. A distinct advantage of this approach is that it allows the pre- 
cise centroiding of the individual galaxies to be included as a free 
parameter in the fitting process, which can compensate for slight 
astrometric differences between the high-resolution and IRAC im- 
ages. However, this method obviously has the disadvantage of rely- 
ing on axisymmetric galaxy models, and is only really suitable for 
deconfusing relatively small areas of IRAC imaging at a time. 

The method adopted in this paper relies on a related, but dif- 
ferent, approach, whereby the actual two-dimensional light distri- 
butions of individual galaxies in the HST imaging are used as the 
model templates. In this approach, SEXTRACTOR is used to pro- 
duce a normalized template of each individual galaxy based on ei- 
ther the J125 or #160 WFC3/IR imaging, and is transformed to 
IRAC spatial resolution via convolution with a transfer function: 



x = E 



[Dij - ELi a * M *j. 



(A2) 



where is the (i, j)-th pixel of the IRAC image, <Jij is the cor- 
responding uncertainty and M t * is the (i,j)-th pixel of the fc-th 
galaxy template. The minimum x 2 occurs when: 
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which can be rearranged as: 
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which describes a set of n linear equations, which can be re-cast as 
a matrix equation: 



Aa = b 
where 
A m k = 



(A5) 



2^ „?. 



and b k 



2-*, rr2. 



From equation A5 it is possible to find a, the n-dimensional vec- 
tor of galaxy template amplitudes, using standard matrix inversion 
techniques. Moreover, the variance of the fitted amplitudes are au- 
tomatically provided by the diagonal terms of the inverse matrix: 



(A6) 



We should note that the basic algorithm outlined here is very sim- 
ilar to that employed by the publicly available TFIT (Laidler et al. 
2007) and ConvPhot (De Santis et al. 2007) software packages. 



APPENDIX B: CANDIDATE PHOTOMETRY, SED FITS 
AND POSTAGE-STAMP IMAGES 

In Tables Bl, B2 & B3 we provide the photometry for each mem- 
ber of our final robust sample. In Figures Bl, B2 & B3 we illustrate 
the results of our SED-fitting analysis by showing the best-fitting 
high-redshift galaxy template, the best-fitting low-redshift alterna- 
tive solution and the dependence of \ on photometric redshift. 
This information is included to allow the reader to judge for them- 
selves the robustness of each high-redshift candidate. In Figure B4 
we show the results of fitting the twenty-one object sub-sample de- 
tected at 3.6/im with SED models featuring nebular continuum and 
line emission. Finally, in Figures B5, B6 & B7 we also provide 
3"x 3"postage-stamp images of each high-redshift candidate in 
the z$ 50 , ^105/^098,^125 & #160 filters. 



WpSF *T = IpSF 



(Al) 



where T is the transfer function and Wpsf and Ipsf are the 
WFC3/IR and IRAC PSFs respectively. The fundamental assump- 
tion behind this technique is that it is possible to reproduce the 
observed IRAC image using a linear combination of these galaxy 
templates. The key advantage of this approach is that, provided the 
astrometry match between the WFC3/IR and IRAC images is suf- 
ficiently accurate that the centroid of each galaxy template can be 
held fixed, the amplitude of each template can be uniquely (and 
analytically) determined by x 2 minimisation, as follows: 
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Table Bl. Basic observational properties of the final high-redshift galaxy sample in the HUDF field. Columns one to three list the candidate IDs and coor- 
dinates. The remaining columns list the photometry of each candidate in the F775W, F850LP, F105W, F160W, IRAC1 (3.6/an) and IRAC2 (4.5/im) filters, 
as measured in a 0.6" diameter aperture, along with their corresponding uncertainties. The magnitudes listed here are derived from the actual fluxes used in 
the SED fitting and are not corrected to total, but have been corrected for galactic extinction and the relative aperture losses between ACS and WFC3/IR. All 
candidates are undetected in filters at shorter wavelengths than F775W, and all detections which are significant at less than the 2a level are listed as 2a upper 
limits. It should be noted that because the IRAC photometry is derived via a deconfusion process the uncertainties and 2a limits are highly position dependent. 
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Table B2. Basic observational properties of the final high-redshift galaxy sample in the ERS field. Columns one to three list the candidate IDs and coordinates. 
The remaining columns list the photometry of each candidate in the F775W, F850LP, F098M, F160W, IRAC1 (3.6^tm) and IRAC2 (4.5/im) filters, as measured 
in a 0.6" diameter aperture, along with their corresponding uncertainties. The magnitudes listed here are derived from the actual fluxes used in the SED fitting 
and are not corrected to total, but have been corrected for galactic extinction and the relative aperture losses between ACS and WFC3/IR. All candidates are 
undetected in filters at shorter wavelengths than F775W, and all detections which are significant at less than the 2a level are listed as 2<r upper limits. It should 
be noted that because the IRAC photometry is derived via a deconfusion process the uncertainties and 2a limits are highly position dependent. 



ID RA(J2000) DEC(J2000) F775W F850LP F098M F125W F160W IRAC1 IRAC2 

ERS.7086 03:32:34.75 -27:40:35.1 >28.28 27.25±°;^ 26.781°;^ 26.89^ 

ERS.6066 03:32:07.86 -27:42:17.8 >28.19 27.27l°;* 3 i 26.55±3;J| 26.68 j 

ERS.9869 03:32:15.40 -27:43:28.6 >28.04 27.52±°;g 26.99±°^ 27.35^ 

ERS.8668 03:32:27.96 -27:41:19.0 >28.41 27.63±g|| 27.1ll° ) ;! 4 27.17^ 

ERS.9100 03:32:20.24 -27:43:34.3 >28.45 27.87l°j® 27.161°,;,!,° 27.24^ 

ERS.7225 03:32:36.31 -27:40:15.0 >28.12 >27.88 27.00l°,!| 26.92j 

ERS.6438 03:32:25:28 -27:43:24.2 >27.99 >27.82 26.771°: 13 26 - 76 - 

ERS.6263 03:32:06.83 -27:44:22.2 >28.17 27.421°;*° 26.7ll°,;J 7 5 26.73^ 

ERS.7776 03:32:03.77 -27:44:54.4 >28.29 27.95l°j° 27.12l°,; 2 ° 27.03 j 

ERS.5847 03:32:16.00 -27:43:01.4 >28.52 27.68l° ) ; 2 2 26.601°'^ 26.63^ 

ERS.8987 03:32:16.01 -27:41:59.0 >28.06 >28.16 27.621 ,;* 6 , 27.22^ 

ERS.3679 03:32:22.66 -27:43:00.7 >28.24 27.28lg;^ 26.19±g;g| 25.94^ 

ERS.7412 03:32:09.85 -27:43:24.0 >28.18 >27.87 26.761°," 26.96j 

ERS.6427 03:32:24.09 -27:42:13.9 >28.38 >27.89 27.141°;^ 26.76j 

ERS.8858 03:32:16.19 -27:41:49.8 >28.12 >27.70 27.47±g'.27 27 - 20 - 

ERS.7376 03:32:29.54 -27:42:04.5 >28.07 >28.19 27.17±g;f| 26.95j 

ERS.8176 03:32:23.15 -27:42:04.7 >28.01 >28.25 27.061 , 22 , 27 - 09 - 

ERS.7672 03:32:10.03 -27:45:24.6 >28.06 >28.01 27.07l°, 2 | 27.01 j 

ERS.7475 03:32:32.81 -27:42:38.5 >28.49 >28.28 27.42±g;|| 26.97j 

ERS.7236 03:32:11.51 -27:45:17.1 >28.07 >27.96 27.67±g;|| 26.93j 

ERS.9041 03:32:23.37 -27:43:26.5 >28.39 >27.84 >28.48 27.23^ 

ERS.10288 03:32:35.44 -27:41:32.7 >28.46 >27.97 >28.76 27.40j 

ERS.8584 03:32:02.99 -27:43:51.9 >28.16 >28.41 >28.52 27.16^ 



12 

.11 


26.99l°; 1 1 I 


26.82l°° 


>26.79 


.08 


26.50lg;l? 


24.86l°;2* 


94 99+ ' 24 


16 
.14 


97 ^'J+0-34 
z;.jj_ Q 26 


> 26.45 


>26.39 


16 
.14 


27 "Sf,+ ' 36 
Z/.3D_ 27 


>26.65 


>26.58 


12 
.11 


27.201°;!* 


26.391°;^ 


26 59+ ' 43 


.08 
.08 


26.821°;!° 


24.88l°i* 


25 23+ - 2 * 


:o° 


ZD --"-0.10 


Z^.6J_ 23 


ic 70+O.58 
zj./j_q 3g 


. 14 
.12 


26 95+ - 18 


26.25l°j9 


>26.77 


12 
.11 


27.281°;!° 


26.361°;^ 


>26.98 


.08 
.07 


26.731°;!^ 


z3.D1_q 42 


>25.65 


16 
.14 


01 17+0.19 
z/ - i/ -0.16 


>25.36 


94 «+°- 35 

Zf .DO_ Q 2 g 


08 
.07 


Z3.yi_ 07 


25.07lg;|f 


95 f,7+ - 50 
ZJ.O/_ 34 


.10 
.10 


26.661";!? 


26.891°;^ 


>26.61 


.10 
.09 


26.87±°;°g 


25 67 +0 3 ° 
Z9.o;_ 23 


25 99+ - 26 


14 
.13 


27.401°;!* 


>27.15 


>26.94 


.10 
.09 


27 05+ ' 27 

Z/.UJ_ Q 22 


96 44+ - 69 
zo.ft_ 42 


>26.41 


.13 
.12 


27.50l°; 2 * 


>27.06 


>26.83 


21 
.18 


27.101°;!* 


>26.94 


>26.60 


.11 
.10 


26.941";!, 1 , 


> 26.47 


>26.34 


.14 
.12 


97 9^+0- 13 
Z '- ZJ -0.12 


> 26.40 


>26.19 


10 
.09 


97 Q1 +0-38 
z '- y -0.28 


>27.01 


>26.82 


.12 
.11 


27.471°;!° 


>27.19 


>26.89 


.22 
.18 


27.08l°;!I 


>26.32 


9S 7Q+0-54 



ERS.8496 03:32:29.69 -27:40:49.9 >28.63 >27.81 
ERS.9923 03:32:10.06 -27:45:22.6 >28.00 >27.92 



27.18l°; 2 ! 



26.82Z 



97 14+0-18 



27.367 



27.41 _ 
27.291 



0.25 
0.16 



25 36+ ' 24 95 9I+ - 37 

ZJ.OO_g 2 q ZJ.5>1_q 2g 



>26.95 



>26.71 



© 0000 RAS, MNRAS 000, 000-000 



A robust sample of galaxies at redshifts 6. 0<z<8. 7 23 



Table B3. Basic observational properties of the final high-redshift galaxy sample in the HUDF09-2 field. Columns one to three list the candidate IDs and 
coordinates. The remaining columns list the photometry of each candidate in the F775W, F850LP, F105W, F160W and IRAC1 (3.6^tm) filters, as measured 
in a 0.6"diameter aperture, along with their corresponding uncertainties. The magnitudes listed here are derived from the actual fluxes used in the SED fitting 
and are not corrected to total, but have been corrected for galactic extinction and the relative aperture losses between ACS and WFC3/IR. All candidates are 
undetected in filters at shorter wavelengths than F775W, and all detections which are significant at less than the 2a level are listed as 2<r upper limits. It should 
be noted that because the IRAC photometry is derived via a deconfusion process the uncertainties and 2a limits are highly position dependent. 
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Figure Bl. SED fits for each member of the final HUDF sub-sample. In each plot the solid (blue) line is the best-fitting z > 6 galaxy SED template, and 
the dotted (red) line is the best-fitting alternative low-redshift solution (z < 2.5). All data points which are detected at less than lcr significance are shown 
as lcr upper limits. In each case the inset panel shows x 2 versus redshift, produced after marginalizing over all other free parameters. The horizontal red line 
in the inset is set a (Xmin + an d indicates the threshold for determining whether or not the competing low-redshift solution can be excluded at the 95% 
confidence level. 
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Figure Bl. Continued. 
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Figure B2. SED fits for each member of the final ERS sub-sample. In each plot the solid (blue) line is the best-fitting z > 6 galaxy SED template, and the 
dotted (red) line is the best-fitting alternative low-redshift solution (z < 2.5). All data points which are detected at less than la significance are shown as 
la upper limits. In each case the inset panel shows \ 2 versus redshift, produced after marginalizing over all other free parameters. The horizontal red line 
in the inset is set a (Xmin ant * indicates the threshold for determining whether or not the competing low-redshift solution can be excluded at the 95% 

confidence level. 
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Figure B2. Continued. 
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Figure B3. SED fits for each member of the final HUDF09-2 sub-sample. In each plot the solid (blue) line is the best-fitting z > 6 galaxy SED template, and 
the dotted (red) line is the best-fitting alternative low-redshift solution (z < 2.5). All data points which are detected at less than lcr significance are shown 
as lcr upper limits. In each case the inset panel shows \ 2 versus redshift, produced after marginalizing over all other free parameters. The horizontal red line 
in the inset is set a (Xmin + an d indicates the threshold for determining whether or not the competing low-redshift solution can be excluded at the 95% 
confidence level. 
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Figure B4. SED fits featuring nebular continuum and line emission for each member of the twenty-one object sub-sample detected at 3.6/im. In each plot the 
thick red line is the best-fitting r = 0.4 Gyr stellar population model (with no nebular emission) and the thin blue line is the best-fitting nebular model (see 
Section 5.3 for a full discussion). All data points which are detected at less than la significance are shown as lcr upper limits. 
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Figure B4. Continued. 
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Figure B5. 3" X 3 "postage- stamp images of the members of the final HUDF sub-sample in the zsso, V105, J125 & Hieo filters (left to right). The postage- 
stamps are orientated such that North is top and East is left. 
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Figure B6. 3" X 3"postage-stamp images of the members of the final ERS sub-sample in the 2850,^098, J125 & -H16O filters (left to right). The postage- 
stamps are orientated such that North is top and East is left 
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Figure B6. Continued. 
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Figure B7. 3"x 3"postage-stamp images of the members of the final HUDF09-2 sub-sample in the 2850,^105, J125 &-H16O filters (left to right).The 
postage-stamps are orientated such that North is top and East is left 
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